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ON  THE  FREQUENCY  FILTERING  OF  TRANSIENT  NOISE  SIGNALS 

Prepared  by: 

Edward  C.  Whitman 


ABSTPACT:  A  class  of  transient  random  signals  is  modeled  as  the  product  of  a 
deterministic,  square  integrable  envelope  function  and  a  Gaussian  random  process 
having  a  well-defined  power  spectrum.  The  passage  of  an  ensemble  of  such  random 
signals  through  a  linear  filter  is  studied  with  particular  emphasis  on  the  mean 
and  variance  of  the  total  output  energy.  It  is  found  that  an  important  role  is 
played  in  these  considerations  by  the  covariance  function  between  values  of  the 
energy  density  spectrum  of  a  sample  function  evaluated  at  different  frequency 
arguments.  Accordingly,  the  form  of  this  function  is  derived  and  portrayed  as  a 
surface  lying  above  a  two-dimensional  frequency  plane.  Examples  of  these  spectral 
covariance  surfaces  are  presented  and  discussed  for  both  rectangular  and  decaying 
exponential  pulses  of  both  broad  and  narrow  band  Gaussian  noise,  and  their  general 
characteristics  are  identified.  Finally,  the  problem  of  idealized  narrow  band 
filtering  is  specifically  approached  and  approximate  expressions  derived  for  the 
mean,  variance,  and  normalized  standard  deviation  of  the  output  energy  of  a  narrow 
band  filter  excited  by  rectangular  pulses  of  narrow  band  Gaussian  noise.  The 
approximate  relationships  of  the  filter  bandwidth,  pulse  duration,  and  underlying 
noise  spectrum  are  explored  for  their  effects  on  spectral  resolution  and  statistical 
stability,  leading  to  an  uncertainty  principle.  The  implications  of  these  findings 
for  spectral  analysis  and  monopulse  signal  processing  are  discussed  in  the  light  of 
this  uncertainty  principle  and  the  limitations  it  imposes  on  the  simultaneous 
precision  of  frequency  resolution  and  spectral  amplitude. 
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Chapter  I 
INTRODUCTION 


In  the  predecessor  of  this  report,  reference  (a),  a  study  was  made  of  the 
statistical  characteristics  of  a  class  of  transient,  random  waveforms  appearing 
similar  to  sonar  reverberation  and  EER  returns,  by  assuming  a  relatively  simple 
mathematical  model  for  these  phenomena  and  following  the  statistical  implications 
of  the  model  to  their  logical  conclusion.  The  intent  of  that  study  was  to  prepare 
the  groundwork  for  a  theory  that  would  provide  some  guidelines  for  predicting  the 
result  of  operating  on  such  waveforms  for  purposes  of  signal  processing  and/or 
measurement.  Also  raised  were  certain  philosophical  questions  about  the  extent  to 
which  accurate  predictions  of  random  transient  phenomena  could  be  made  and  about 
the  measurement  interpretation  of  the  harmonic  spectrum  of  such  signals.  This 
report  deals  with  a  considerable  amplification  of  the  latter  point. 

The  model  treated  in  reference  (a)  is  shown  in  Figure  1.  It  consists  of 
generating  an  ensemble  of  random  transient  waveforms  by  the  simple  expedient  of 
multiplying  a  continuing,  stationary,  zero-mean  random  process  by  a  deterministic, 
square-integrable  transient  which  is  taken  to  be  zero  for  t  <  0.  By  considering 
an  infinite  ensemble  of  stationary  random  processes  and  a  single  transient  "envelope 
waveform,"  one  generates  an  infinite  ensemble  of  non-stationary  random  transients 
each  differing  from  the  other  in  detail,  but  all  characterized  by  the  same  time- 
varying  variance,  as  the  following  shows. 

If  we  denote  a  typical  output  ensemble  member  by  s(t),  the  continuing  station¬ 
ary  input  process  by  n(t),  and  the  multiplying  transient  waveform  by  e(t),  then 
evidently 


s ( t )  =  e(t)  n(t) 


(1) 


If  a2  is  the  variance  of  n(t),  and  o2  that  of  s(t),  then  since  e(t)  is 
assumed  deterministic  and  hence  does  not  enter  into  the  averaging  process  for  the 
variance, 


o2(t)  =  e2(t)a2  (2) 

s  n 

which,  as  noted  above,  is  a  function  of  time. 

For  each  member  s(t)  of  the  ensemble  of  random  transients,  one  can  compute 
the  time  autocorrelation  function  defined  by 


FIG .  1  A  MODEL  FOR  THE  GENERATION  OF  RANDOM  TRANSIENTS . 


00 


(3) 


4>  (r)  =  /  s(t)  s(t+T)dt 

SS  0 


and  the  Fourier  integral 


/  s  (t)  e-^  dt 
o 


(4) 


Thd  Energy  density  spectrum  of  the  ensemble  member  is  then  given  by 


#ss(u)  =  2n  S(u)  SW  =  211  |S(m)  | 2 


(5) 


where  the  superscript  bar  denotes  complex  conjugation.  From  the  Wiener-Khintchine 
relation,  we  know  also  that 


_1 

2n 


00 


/ 

00 


rss 


(t)  e  dr 


(6) 


- 


which  is  to  say  that  the  time  autocorrelation  function  of  equation  (3)  and  the 
energy  density  spectrum  of  equation  (5)  constitute  a  Fourier  transform  pair.*  If 
the  time  autocorrelation  function  and  energy  density  spectrum  are  computed  for  each 
of  a  large  number  of  sample  waveforms,  they  will  turn  out  to  be  different  from 
sample  to  sample  since  they  result  from  operations  on  randomly  different  waveforms  • 
Thus,  for  every  value  of  their  arguments,  these  quantities  are  random  variables 
which  may  be  considered  statistics  of  the  ensemble.  It  is  not  possible  to  speak  of 
a  single  deterministic,  predictable  autocorrelation  function  or  spectrum  valid  for 
the  entire  ensemble,  since  each  ensemble  member  will  yield  a  different  function 
which  can  not  be  predicted  with  exactitude  before  the  event. 


Since  the  time  autocorrelation  function  and  energy  density  spectrum  of  a  random 
transient  are  random  variables,  parameterized  in  a  sense  by  their  respective  argu¬ 
ments,  it  becomes  interesting  to  study  the  expected  values  of  these  variables  and 
also  their  variances.  In  practical  terms,  the  expected  values  represent  the  average 


*NOTE ;  The  convention  of  reference  (a)  will  be  followed  here  in  that  autocorrelation 
functions  and  spectra  for  "energy"  (i.e.,  transient)  signals  will  be  denoted  <j>  and 
$  respectively  (using  upper  and  lower  case  phi) ;  whereas  autocorrelation  functions 
and  spectra  for  "power"  (i.e.,  continuing)  signals  will  be  respectively  denoted  \p 
and  ¥  (upper  and  lower  case  psi) . 
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at  a  point  obtained  over  a  large  number  of  samples,  and  the  variances  indicate  the 
extent  to  which  the  spectrum  and  autocorrelation  of  a  single  sample  can  be  expected 
to  depart  from  the  average  at  a  point.  In  reference  (a),  these  means  and  variances 
were  derived  both  in  general  and  in  particular  simple  cases  where  the  spectral 
character  of  the  continuing  process  and  the  form  of  the  envelope  transient  were 
specified.  The  principal  results  were  the  following: 

If  s(t)  =  e(t)n(t)  and  <}>ss(t)  and  $ss(m)  are  given  by  equations  (3)  - 
(6) ,  then 


E[<j>  (x)]  =  (t)<6  (t)  (7) 


<L  (x)e'jt0X  dx  (8) 

06 


UU 


and 


E[$>ss(a>)] 


2n  _  f  ^nn^ 


where  E[ * ]  denotes  an  ensemble  average,  $nn(x)  th6  autocorrelation  function 
of  the  continuing  stationary  process  and  d>ee Ct)  is  the  autocorrelation  of  the 
envelope  waveform.  The  second  of  these  relations  can  be  expressed  as 


00 

E[$ss(oj)]  =  /  ^(w’)  $ee(w  ~  w’)d(i)’ 

—  00 


-  V(w)  ®  *ee<“> 


(9) 


where  Tnn(w)  is  the  power  spectral  density  of  the  noise  process,  ^eR(w)  is  the 

energy  density  spectrum  of  the  transient  envelope,  and  represents  the  operation 
of  convolution.  The  expected  value  of  <j>ss(x)  is  thus  seen  to  be  the  product  of 
the  autocorrelation  functions  of  the  two  signals  whose  product  forms  the  random 
transient,  and  not  surprisingly  the  expected  value  of  the  energy  density  spectrum 
turns  out  to  be  given  by  the  convolution  of  the  power  density  spectrum  of  the 
continuing  stationary  process  and  the  energy  density  spectrum  of  the  envelope 
waveform. 

If  n(t)  is  further  restricted  to  be  a  Gaussian  random  process,  the  following 

expressions  for  the  variances  can  be  derived: 

00 

Var [<f>ss (t) ]  =2  /  <j)pp(u,  x)  [^n(u)  +  ^nn^u+T)^nn^u-T^ ^  du  (10) 


where  d>pp  (vi ,  x) 


is  the  autocorrelation  function  of 


P(t,  x) 


e(t)  e(t-x) 


(11) 
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that  is 


4>pp(u,  t) 


p(t,  t)  p(t  +  u,r  )  dt 


Var[4>sg(w)] 


=  {E[$ss(u>)]}2 


w  w 

f  f  e(t)e(u)^  (t-u)e  ^u^t+u^dt  du|2 


The  study  reported  here  is  an  extension  of  the  approach  described  in  reference 
(a)  and  concentrates  on  the  energy  density  spectrum  of  random  transient  signals. 

In  particular,  it  will  be  aimed  at  predicting  the  result  of  frequency  filtering 
such  waveforms  and  at  the  somewhat  broader  question  of  describing  the  distribution 
of  energy  over  the  spectral  band  of  the  signals.  If  one  considers  a  narrow  band 
filter,  for  example,  and  passes  an  ensemble  member  through  it,  the  output  energy 
observed  is  an  indication  of  the  "strength"  of  the  signal  in  the  particular 
frequency  band  passed  by  the  filter.  Upon  repeating  the  experiment  with  other 
ensemble  members,  however,  it  will  be  found  that  in  general  this  output  energy  will 
vary  from  sample  to  sample  and  is  thus  a  random  variable  parameterized  not  only  by 
the  characteristics  of  the  random  transient  ensemble,  but  also  by  the  bandwidth  and 
center  frequency  of  the  filter.  It  is  of  interest  to  know  the  expected  value  and 
variance  of  this  output  energy  as  an  aid  to  measurement  planning  or  performance 
prediction,  and  this  is  the  task  undertaken  in  the  following  pages.  It  will  be  seen 
that  the  resulting  analysis  bears  a  strong  resemblance  to  that  underlying  traditional 
power  spectrum  measurement,  which  indeed  can  be  viewed  as  a  special  case  of  the  more 
general  problem  treated  here.  One  of  the  goals  of  the  present  study  is  the  develop¬ 
ment  of  a  proper  interpretation  of  the  energy  density  spectrum  of  a  transient  signal 
and  its  relationship  to  physical  measurements  of  the  signal  energy. 
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Chapter  II 

FREQUENCY  FILTERING  OF  RANDOM  TRANSIENT  WAVEFORMS 


To  study  the  effects  of  frequency  filtering  on  random  transient  signals,  one 
need  only  make  a  fairly  straightforward  application  of  linear  system  theory.  If 
the  filter  of  interest  is  modeled  as  a  linear  system  with  impulse  response  h(t), 
it  is  well  known  (reference  (b) ,  Chapter  13)  that  if  s(t)  is  the  input  to  the 
filter,  then  the  output  f(t)  is  given  by 


f(t)  **  /  s(x)  h(t-x)  dx  (14) 

o 

which  is  the  convolution  of  the  input  signal  and  the  impulse  response  of  the  filter. 
In  terms  of  energy  density  spectra  this  becomes 

$ff(w)  =  |H(a>)  | 2  $ss(w)  (15) 


where  $ss(w)  and  £  (to)  are  the  input  and  output  spectra,  respectively,  and 
H(oj)  is  the  system  function  of  the  filter,  given  by  the  Fourier  transform  of  h(t): 

CO 

H(w)  =  /  h(t)  e'^t  dt  (16) 

—  CO 

Since  $ss(w)  is  a  random  variable,  it  should  come  as  no  surprise  that  for 
every  value  of  to,  ^^(w)  is  a  random  variable  also.  From  equation  (15), 
however,  it  is  immediately  obvious  that 


E[4>ff(w)]  =  |H(o>)  [ 2  E[4'sg (oj)  ] 


=  |H(w)  I2  [^(to)  ®  ^(w)] 


(17) 
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where  E[S>  (to)]  is  given  by  equations  (8)  or  (9).  Similarly, 

5b  , 

Var  [$>ff(a))]  =  | H(w)  | 4  Var[$sg(u)]  (18)  ; 

with  Var  [ _  (to)  ]  given  by  equation  (13).  Since  the  second  term  of  equation  (13)  ; 

SS 

is  always  positive,  it  is  apparent  that  ' 

Var  [4>ss(u)]  >  (E[<J>gs(u))]}2  (19) 

and  in  turn  that  i 

Var  [$ff(a>)]  >  (E[$ff(o))]}2  (20)  j 

This  is  to  say  that  the  standard  deviation  of  $££((o)  will  always  be  larger  than 
the  mean  of  $££(w),  regardless  of  H(to),  which  implies  that  the  distribution  of 
4>££ (co)  is  rather  broad  for  every  o).  In  practical  terms,  this  means  that  in 
calculating  ^^(w)  for  a  number  of  sample  functions  and  comparing  the  results  for 

some  given  w,  a  large  variability  will  be  found.  If  one  is  seeking  a  statistic  of 
the  random  transient  ensemble  that  remains  fairly  constant  from  sample  to  sample, 
the  value  of  $ff(a))  at  a  point  clearly  falls  far  short  of  the  ideal.  The  situa¬ 
tion  here  is  rather  similar  to  that  faced  in  attempting  a  spectral  analysis  on  the 
basis  of  the  so-called  "periodogram”  (reference  (c) ,  pageo  107-108)  which  can  be 
shown  to  have  little  statistical  reliability  as  an  estimate  of  the  deterministic 
spectrum  of  a  stationary  continuing  stochastic  process,  unless  a  number  are  computed 
and  ensemble  averaged. 

To  gain  a  possibly  more  reliable  statistic  for  the  ensemble  of  random  transients  : 
processed  through  the  filter  H(w) ,  we  will  consider  ’’ntegratong  the  c.tput  spectrum 
$££  (o>)  over  all  to  to  obtain  an  expression  for  t\a  tofal  erargy  in  f^t).  An 
important  special  case  will  be  that  in  which  |H(w) |  is  very  narrow  compared  to  the  ? 
spectral  width  of  s(t),  in  which  event  the  filter,  with  a  tunable  center  frequency,  ! 
provides  in  some  sense  an  empirical  measurement  of  4>ss(oj).  It  will  be  important  to  j 
examine  the  statistical  reliability  of  such  measurements,  i.e.,  the  extent  to  which  \ 
a  single  measurement  accurately  reflects  the  ensemble  as  a  whole.  This  evidently  1 

has  an  important  bearing  on  the  validity  of  attempting  to  measure  or  predict  the  ’ 

spectrum  of  a  sample  random  transient  to  arbitrary  precision. 

A  diagramatic  representation  of  the  problem  of  interest  is  shown  in  Figure  2.  • 

Let  F(H)  represent  the  total  energy  of  s(t)  remaining  after  it  is  passed  through 
the  filter  H(w).  In  terms  of  f(t),  the  output  signal, 

OO 

F(H)  -  /  f2(t)  dt  (21) 


7 
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FIG. 2  LINEAR  FILTERING  OF  RANDOM  TRANSIENT  WAVEFORMS. 
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and  in  terms  of  the  various  spectra, 

oo  oo 

F(H)  =  /  $££(u))do>  =  /  |H(<o)  |2  $ss(a))dti) 

—  CO  —00 


(22) 


F(H)  is  a  random  variable  whose  mean  is  given  by 

00 

E[F(H)]  =  /  |H(w)|2  E[$ss((i>)]dio 

—00 


(23) 

00 

=  /  ] H(o>)  ! 2  [*^(<*0  OD  4>ee(oj)  J d(0  by  equation  (9) 

—  CO 


To  find  the  variance  of  F(H),  we  will  use  the  relationship 


Var[F(H) ]  =  E[F2(H)]  -  {E[F(H)]}2 


(24) 


CO  00 

F2(H)  =  I  | H (to)  | 2  4>ss(o))do)  f  |h((jj’)  I2  <J>SS (w* )da)' 

-00  — oo 


(25) 

CO  00 

=  f  f  |H(u)  I2  |H(u)’)|2  <J>ss(to)  4>ss (a) 1  )du>da) ' 

—  CO  —00 


Taking  the  expected  value  yields 

CO  00 

E[F2(H)]  =  S  /  |H(w)  ] 2  |H(w')|2  E[4>ss(oj)  $ss(o)')]da)da)' 

—  00  — oo 


Now  using  equations  (23)  and  (24) , 


00  00 

Var[F(K)  ]  /  /  |H(u>)|2  |H(m')|2  E[$ss(co)  4>ss (u»f  )]d<od(o' 

—00  —00 


(26) 


CO  00 

-  /  |H((0)|2  E[<J>ss(o))]do)  /  |H(o)')|2  E[$ss(w')]dw' 


(27) 
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A  little  manipulation  then  gives 

00  00 

Var[F(H)]  -  /  /  |H(o>)|2  |H(o)t)|2  {E[$8s(u)$88(ior)] 

—00  —oo 


-  E[4>bs(w)-]  E[$ss ] >  dudio’  (28) 

The  quantity  in  the  brackets  can  immediately  be  identified  as  the  covariance 

of  $  (to)  and  (wr)  since  it  can  be  easily  shown  that 
88  88 


Cov(x,  y)  »  E(xy)  -  E(x)E(y)  (29) 

(see  reference  (d),  Ch.  8,  Sec.  2).  We  may  thus  finally  write  that 

00  00 

Var[F(H)]  =  f  f  |h(o))|2  |H(o>f )  |2  Cov[$ss(w),  ^(w^ldudu' .  (30) 

—  00  —00 


At  this  point,  it  is  well  to  pause  for  an  interpretation  of  the  results 
attained  thus  far.  In  the  situation  considered  here,  sample  members  of  an  ensemble 
of  random  transients,  generated  according  to  the  description  in  Chapter  I,  are 
passed  through  a  linear  system  represented  by  H(w).  As  noted  above,  each  sample 
input  waveform  will  have  a  different  computed  spectrum  $ss(w),  and  hence  each  output 
waveform  will  have  a  correspondingly  different  spectrum  ^(w).  The  quantity  of 

interest  here  is  the  total  energy  of  the  filter  output,  denoted  F(H)  and  given  by 
the  integral  of  (u)  over  all  u>.  F(H)  is  a  random  variable  whose  mean  is  given 

by  equation  (23) .  This  is  the  result  to  be  expected  in  averaging  the  total  output 
energy  over  a  large  number  of  sample  transients.  Graphically,  E[F(H)]  can  be 
visualized  as  the  area  under  the  curve  formed  by  plotting  E[$ss(u))]  over  all 

frequencies  and  then  weighing  by  the  function  |H(to)  | 2 .  The  square  root  of  the 
variance  of  F(H) ,  namely  the  standard  deviation  ,  is  an  indication  of  the 

spread  uC  the  probability  distribution  of  the  total  output  energy  about  its  mean. 
o2p(H)  is  given  by  equation  (30)  above  and  can  be  interpreted  as  the  volume  under 

a  surface  defined  over  the  w  -  u*  frequency  plane.  The  height  of  this  surface 
above  the  plane  is  given  by  the  product  of  Cov[4>  (w) ,  3>  (co1)]  and 

| H(tu)  | 2  )h((o')|2.  Alternatively,  one  can  visualize  a  "spectral  covariance  surface" 
whose  height  above  the  w  -  u’  plane  is  given  by  the  function  Cov[$>ss(io) ,  $ss  (to T )  ] . 
The  variance  integral  can  be  considered  to  yield  the  volume  under  such  a  surface 
after  each  volume  element  has  been  weighted  by  |h(w)|2  | H (to T )  | 2 .  The  spectral 
covariance  surface  thus  plays  a  central  role  in  the  study  of  the  frequency  filtering 
of  an  ensemble,  particularly  when  |h(oj)  |  is  relatively  constant,  and  its  properties 
will  be  examined  in  the  next  chapter. 
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Chapter  III 

SPECTRAL  COVARIANCE  SURFACES  FOR  RANDOM  TRANSIENT  ENSEMBLES 


In  this  chapter,  we  shall  study  the  properties  of  the  covariance  of  (or 
between)  values  of  $  (a) )  for  two  arguments  and  M  .  From  elementary 

statistics, 

Cov[$ss(u),  $ss(u’)]  »  Ef$CQ(o))  -  E[$ss(to)]  Et^Cco')] 


‘ss"~'  ss 


If  the  Fourier  Transform  of  s(t)  is  given  by 


S(u)  =  -i—  *  s(t)e~jUtdt, 

2n  „ 


*sg(w)  =  2n  s(»)  s(w)  -  2n|s(u)|2 


Recalling  that  s (t)  a  e(t)n(t), 


$ss(u)  =  ~~zT  J  e(t)n(t)e"ja)t  dt  /  e(x)n(x)ej“X  dx 


00  00 

~  f  f  e(t)e(x)n(t)n(x)e  t  x^dt 


Changing  one  of  the  variables  of  integration  to  a  ■  t-x  yields 

00  00 

$  (u)  -  ^  f  !  e(t)e(t-a)n(t)n(t-a)e  &  dt  da 

38  2n  -00  0 


In  taking  the  expected  value,  we  use  the  fact  that  n(t)  is  stationary  such  that 


E[n(t)n(t-a)  ]  *  ^(a) 


11 
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and  that  /  e(t)e(t-ra)dt  5  <j>  (a) 

ee 


Thus,  »  co 

El$gsCu)]  -  2n~  ^  e(t)e(t~a)$nn(a)e  dt  da 

**oo  0 


00  ♦ 

_i_  /  (a.)  <a)e  Wa  da 

2n  -a,  ee  nn 


which  was  previously  given  ex  cathedra  as  equation  (8). 

Similarly,  the  product  <I>SS Cto)  $>ss(<i)’)  may  be  written  as 


*ss *ss(a),) 


CO  CO  CO  CO 

f  f  f  f  e(t)e(x)e(u)e(v)n(t)n(x)n(u)n(v)e  X  e-^  ^u-v^dt  dx  du  dv  (39) 


oooo 


If  n(t)  is  a  Gaussian  random  process,  it  may  be  shown  (reference  (e) ,  page  93) 
that 

E[n(t)n(x)n(u)n(v )]  *  f  (t-x)f  (u-v) 

nn  rnn 


+  ^nn(t-uHnJX~v)  +  ^nn(t~v^nn(x-u) 


and  therefore, 

^ssfo^ssfo')] 


i  /  /  /7e(t),M,(*(,)t  at  dx  du 

r\  r\ 


oooo 


00  00  00  00 

+  IP-  f  f  f  1  e(t)e(x)e(u)e(v)^nn(t-u)i|)nn(x-v)e~^(t"x)e”:5a),(u"v) 


oooo 


dt  dx  du  dv 


w  w  w  oo 

f  TjfZ  f  1  1  !  e ( t) e (x) e  (u)  e  (v) i/>  ( t-v) ij>n  (x-u) e~^ 03  ( t_x) e-3 u '  (u-v)  dt  dx  du  dv 
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which  will  be  abbreviated  as 


E[*SSM  *„(«')]  -  Ij.  +  I2  +  x3  (42) 

where  Ip  I2,  and  I3  represent  the  three  integrals  of  equation  (41). 
can  immediately  be  written  as 

00  00  .  .  00  00 

X1  *  7n7  J‘  *  ettJeCxJ^ft-xJe  U  t  X)dtdx  I  /  'u^v^^u-NOe-^' (U-V) 


0  0 


0  0 


dudv  (43) 


Considering  the  expected  value  of  both  sides  of  equation  (34)  shows  that 


h  =  E[*8S(a))]  E[$ss(a),)] 


Continuing 


(44) 


(45) 


T  .  1  "  "  ,  -j  (Wfrho'u)  “  “  ,  w  ,,  ,  ,  j  ((OX+Oj’v)  . 

I2  4ii7  f  f  e(t)e(u)4»nn(t-u)e  dtdu  /  /  e(x)e(v)ijinn(x-v)e  dxdv 


0  0 


0  0 


Noting  that  the  two  double  integrals  are  complex  conjugates  of  each  other,  one 
writes 

00  00 

,  J  - - 'a+aJ 

rnn' 


-Iff  e(t)e(u)i|)T1n(t-u)e"^wt+t°  u^dtdu 
211  0  0 


(46) 


Now  using  the  change  of  variables  a  =  t-u, 


X2 


f  f  e(t)e(t-a)ip  (a)e"'ja>t.  e_JW  (t_a)  dadt 
211 _  nn 


00  0 


(47) 


2H  /  9  (  /  e(t)e(t-a)e  3  ^ 01+10  dt}da 


A  further  simplification  can  be  wrought  by  using  the  function  defined  in 
equation  (11) : 


p(t,  x) 


e(t)e(t-x) 
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This  has  £  Fourier  transform  :  '>• 

j,r-*  -  ->  -  ia»  •  -  *  *  *.  '  ,•  l  ' :  ' 

.  ?(«>,  x)  “  5§“  J"  pCt,  a£)e”^Ut  dt  ,£.  (48) 

^  .  ;0  -  - 

12  may  now  bo  expressed  as  '  • 

l£  *  J^-  ^  T(>j)  +  ttii}e2u*x  dx|2  (49) 

.  .  _  A  IV  "  ,  *  ’  "  ,  ‘  ' 

Similarly  it  may  be  shown  that 

l  *”  I2 

I-  -  -A-  s  ^(x)  P(u  -  u\  x)e  2  dx]  <5D) 

1 2n  — la 

Note  that  It  As  merely  equal  to  I2  with  w*  replaced  by  (-w*)..  Finally,  by  ’v 
equations  (31)  and  (44),  • 

Cov[$£8(u;),  «S8(u')3  -  El*e,»  $gB(m*)]  -  EIObb(o))]  Flf6B(o)’)] 


-  I2  +  12  *  I3  -  IA 

-  h  +  h 

Cov  I4'6S(w).  $SE(w’)J  - 

to 

I  W~-J  P'“  +  “*♦  x)e3“  X  dxj 

-  •  1  *  2 
+  A-  /  It'jjjjOO  ?(<*>  T  (D1*  x)e  x  dx! 


<51) 


(52) 


where  P(uj,  x)  is  given  by  equation  (48). 


Some  fundamental  properties  of  the  spectral  covariance  function  can  be 
appreciated  at  once  by  considering  the  form  of  equation  (52).  First  we  note  that 
due  to  the  very  definition  of  covariance,  the  function  will  not  be  changed  In 
value  if  w  and  w ’  are  interchanged.  From  equation  (52),  it  is  apparent  that 
replacing  to'  by  -<o’  will  similarly  have  no  effect  on  the  value  of  the  function. 
From  these  two  facts,  it  becomes  obvious  that  u  and  -u  are  also  interchangeable 
and  thus  that  Cov  [$gB(w),  '5>ss(w,)3  is  an  even  function  of  u>  and  u’  separately 

and  also  a  symmetrical  function  of  u>  and  u’ .  Geometrically,  this  means  that 
the  covariance  function,  defined  as  a  surface  over  the  u  ~  w’  plane,  must  be 
symmetrical  about  the  vertical  planes  defined  by  the  equations  ui  »  u>’  and 
w  *  -u* ,  and  preserve  the  same  shape  (e-xcept,  of  course,  for  a  rotation)  above  aJ. 
four  quadrants  of  the  w  -  m'  plane.  Cov  [C’BS(«),  $  (u>')3  can  thus  be  said  to 
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possess  a  cruciform  symmetry  about  vertical  planes  inclined  45°  with  respect  to 
those  defined  by  the  equations  m'  =  0  and  m  =  0.  This  is  shown  in  Figure  3. 
Evidently,  it  is  necessary  to  study  the  function  only  in  a  single  quadrant  to 
visualize  it  everywhere.  Another  interesting  property  of  Cov  [$c_(aj),  $  (to1)] 

that  emerges  directly  from  equation  (52)  is  that  it  is  always  positive.  This  was 
by  no  means  obvious  at  the  outset,  since  two  random  variables  may  in  general 
exhibit  negative  covariance  even  if  it  can  be  shown  that  they  are  separately 
always  positive  (as  $  (io)and  (co’)  are  here).  Values  of  a  sample  spectrum  for 

SS  ss 

two  different  arguments  are  therefore  always  positively  correlated.  It  is  not 
really  clear  at  this  point  what  deeper  significance  this  may  hold. 

To  continue  this  investigation  of  the  form  of  the  spectra]  covariance  surface, 
it  is  necessary  to  turn  to  specific  expressions  for  the  envelope  function.  Accord¬ 
ingly,  the  two  forms  treated  in  reference  (a)  will  be  considered:  a  rectangular 
pulse  of  duration  T  seconds,  and  a  decaying  exponential  pulse  having  a  time 
constant  of  T  seconds. 


For  the  first  of  these,  shown  in  Figure  4, 


e(t)  =1,  0<  t<T 

=  0,  elsewhere 

From  equation  (11),  for  0  _<  x  _<  T, 

p(t,  x)  =  1,  x<  t  <T 

=  0,  elsewhere 


and  for  -T  <  x  <  0 


p(t,  x)  =1,  0jf_t<_T  +  x 


=  0,  elsewhere 


(53) 


(54) 


(55) 


Evidently  for  |x|  >  T,  p(t,  x)  =  0  everywhere.  The  function  p(t,  x)  is 
plotted  in  Figure  4  as  a  function  of  t  with  x  treated  as  a  parameter.  Now 
taking  the  Fourier  transform  with  respect  to  t  as  in  equation  (48)  yields 


P j^(w ,  x) 


(56) 


for  -T  _<  x  £  T  and  -»  <  w  <  ». 

where  the  subscript  R  on  P_ (u> ,  x)  denotes  the  rectangular  pulse  case. 
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When  |x|  >  T,  PR(<o,  x)  =  0.  for  all  w. 


For  the  exponential  envelope  case,  shown  in  Figure  5, 


e(t)  -  e  t/fT,  t  >  0 


t  <  0 


Accordingly,  for  x  0, 


and  for  x  <  0, 


\  x/T  -2t/T 
p(t,  x)  =  e  e  ,  t  >_  x 


,  t  <  x 


p(t,  x)  =  ex/T  e"2t/T,  t  >_  0 


*  0  ,  t  <  0 

From  equation  (48),  with  Pg(w,  x)  denoting  the  Fourier  transform  appropriate  to 
the  decaying  exponential  case, 

-x/T 

P  x')  „  1  e _  -jwx 

’  }  2n  2  e  ,  -»  <  a)  <  »  (60) 

~  + 

for  x  >  0 


PE(w>  x) 


JL _ 

2F  2 


— 00  <  Ot)  <  °° 


+  jlO 


for  x  <  0 


Having  derived  expressions  for  P(w,  x)  in  the  two  envelope  cases  of  Interest, 
it  is  now  possible  to  continue  the  development  by  substituting  these  into  equations 
(49)  and  (50)  for  I^Ga,  cor)  and  u>').  Again  using  the  subscript  R  to  denote 
the  rectangular  envelope  case, 


Iob(w,u)’ )  -  /  tp 


ooizlsl 

nn1'  J  211  ,  .  i \ /’T—  I x  I  \  e  3  \2  /  e 


(*t,')GdEL)  e 


-j (w+w 1 ) 


3-  f\  w  ids! 

-t  nn  2n  (u+u')(EET) 


-j(tO-O)’)  X 


18 
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Noting  that  the  integrand,  except  for  the  exponential,  is  an  even  function  of  x, 
and  using  the  Euler  relations,  one  writes 

f  / T-x\  2 

I2R(aJ,aj’)  =*  e~3(wfw,)-|-  /  ^nn(x)  (T-x)~~~~~~ — cos  (ara>,)-J|-  dx  (63) 

0  (***>')  ("TV 


_  *  i  s  T 

The  magnitude  of  e  J  v  J  2  is  unity,  and  the  integral  is  real,  so 


T  sinCoj+a)’)  (  T~x  ) 

T  f  ,  v  1  /t|»  (x)  (T-x) -  \  ?  .  I  COS  (o)-U  * )— r-  dx 

V“’“ )=  -w  [o  nn  2  J 


Recalling  that  13(a),  o>')  =  I2(m,  -a>'),  we  have 


(¥) 


Io„(a),uT)  =  sinCowD*)  l  o  y 

3R  JT2^  /  ifinn(x)  (T-x) - - - 7-3  cos(uH-w,)^2—  dx 

(w-w*)  j 


The  sum  of  equations  (64)  and  (65)  gives  the  covariance  function  for  the  rectangular 

envelope  case  when  i|)  (x)  is  the  autocorrelation  function  associated  with  the 

.  ,  nn 

input  noise  process. 

For  the  exponential  envelope  case,  the  same  procedures  yield  the  following 
results : 


(-f-)2  +  (uHV) 


2  /  ^nn(x)e 

1  o 


-x/T[e-jU'x  +  e-jc 


(66a) 


(-|-)2  +  (c^')? 


/  ip  (x)e~x/T[ejw'*  +  e"jl0X]dx 


(66b) 


or  defining 


I  =  /  ib  (x)e  [cos  u'x  +  cos  rnx)]  dx 

a  nn 

o 


(67a) 
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(67b) 


(67c) 


(68) 


Further  progress  toward  a  description  of  the  spectral  covariance  surfaces  un¬ 
fortunately  requires  the  use  of  specific  expressions  for  the  input  noise  autocorre¬ 
lation  functions.  This,  to  be  sure,  destroys  the  generality  of  the  results  (which 
to  some  extent  has  been  sacrificed  already  by  using  specific  envelope  functions), 
but  by  carefully  choosing  the  noise  spectra  of  interest,  one  can  arrive  at  fairly 
general  conclusions  about  the  form  of  the  covariance  surface  for  an  arbitrary  noise 
type.  In  reference  (a),  two  examples  of  idealized  noise  spectra  were  treated,  and 
the  same  will  be  used  here.  These  are  the  rectangular  broad  band  snectrvtm,  based 
on  zero  frequency,  and  a  rectangular  narrow  band  spectrum  centered  at  some  coo. 

The  two  spectra  are  portrayed  in  Figure  6. 

The  first  can  be  described  mathematically  as 


"  -x/T 

lb  5  /  ^nn(x)  e  [sin  co’x  +  sin  cox]  dx 


00  -x/T 

I  =  f  ^nn(x)  e  [sin  <o!x  -  sin  cox]  dx 
o 


covE[$sS(co),  *S8(<o')]  -  4 p-  JJf 


(j)  +  (tu  +  (0*) 


[I  2  +  I*2] 


m2  1  a  b 


+  &  (Tf+  ZT-  .•)*  [I°2  +  *'2] 


¥nn(w>  =  No>  -a)2  -  “  -  “2 
=  0,  elsewhere 

The  corresponding  autocorrelation  function  is 


(69) 


sin  co2t 

\b  (t  )  =  2N  co«  -  ,  -oo  <  t 

rnn  o  2  co  t 

2 


<  00 


(70) 


The  second  spectrum  is  given  by 


«nn(w)  =  Nq,  o'j  <  |coj  <.  co2 


(71) 


0,  elsewhere 


NOLTF  68-124 


.  „  .  sin (Am) t 

ipnn(r)  =  4AcoN0  - (^7-  cos  woT  ’  <  T  <  “ 


where 


Aco  = 


6)  2  — 


,  the  half  bandwidth 


w,  to 
1  +  2 


,  the  center  frequency 


(72) 

(73a) 


(73b) 


These  expressions  have  been  used  in  equations  (64) ,  (65) ,  (67) ,  and  (68)  to 
derive  integral  representations  for  the  spectral  covariance  surfaces  of  both 
rectangular  and  decaying  exponential  bursts  of  both  broad  and  narrow  band  noise. 
These  expressions  have  in  turn  been  evaluated  on  the  IBM-7090  computer  at  NOL  to 
gain  an  understanding  of  the  form  of  the  surface.  Indeed,  much  of  the  manipulation 
leading  up  to  the  final  expressions  for  Cov [ <J>sg (m) ,  $ss(u)T)]  was  specifically 

intended  to  reduce  the  problem  to  one  of  evaluating  Gaussian  quadratures  -  a  task 
amenable  to  a  computer  algorithm.  We  will  now  consider  several  of  the  cases  in 
detail. 


A.  Rectangular  Pulse  of  Low-Pass  Broad  Band  Noise 

For  the  rectangular  function  described  in  equation  (53)  it  is  easily 
shown  that 

<j>  (x)  =  T  -  |  t  | ,  -T  <  r  <  I  (74) 

ee 

=  0,  elsewhere 


Thus,  from  equation  (38), 


E[*ss(o>)] 


1_ 

2n 


f  *eeU)  *nn(x)e 

>00 


-jwx 


dx 


(75) 


~  j 1  /  (T-x)  \p  (x)  cos  wx  dx 

0 

and  from  equations  (64)  and  (65), 
c°v[ 4>ss (to)  ,  ♦„(«')]  = 


1 


r  T  sin(u>+u)')  (Xp) 

/  P _ (x)  (T-x) 


nn 


(urt-u)' 


— r=— ~  cos  (w  -  u>’)  9  dx 


sin(o) 


!  *nn(x)  (T-x)  cos  (»-»')  ^  dx 


(76) 
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It  Is  also  convenient  to  recall  that 


Var[$  (to)]  =  Cov($  (10) ,  $  (w)] 

ss  ss  ss 


or  from  equation  (76)  with  to  =  to', 


Var[4>  ..(to)]  = 

b  5 


It2 


sin  to(T-x) 


1  .WHO 


dx 


f  ib  (x)  (T-x)  cos  tox  dx 
nn 


(7 


(7 


■  I_  (w,  to)  +  {E[$  (to)]}2 

zK  SS 

where  in  the  last  expression  we  have  used  equations  (64)  and  (75). 

To  obtain  the  final  expressions,  we  substitute  ^nn(x)  from  equation  (70) 
above  into  equations  (75) ,  (76) ,  and  (78) , 


EI*ss<“>J 


n  Mi-$) 


a^x 


cos  wx  dx 


(7 


Cov[$gs(w),  *gg(to*)J 


4N  2w22T2 

- x 


T  sin  wox 

Ml  -  i) 


to2x 


sin(to-to')  (^2^-) 

(N-a.'jfe)  <“  +  M) 


dx 


7- 


< 


f  a  _  x^  sin  “2X 

T  " 


sin(to+a)')(^—^) 


to  x 


( .  .  i  \  (  T~x  N 


cos  (to  -  to')(|jdx 


Uc 
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Var[$ss(o))]  =  {E[$ss(a))]>2 


(81) 


4N026J22T2 

n2 


T 

S  tt-s 

o  i 


sina^x  sinco(T-x) 
(i)2x  o)(T-x) 


2 


For  greater  generality  of  presentation,  these  expressions  will  be  normalized 
and  parameterized  by  defining  the  following  quantities,  as  in  reference  (a): 


MpT 


q  =  — 2- 

4  2n  > 

the  number  of  periods  of  the  highest  noise  frequency 
contained  in  a  pulse  length. 

(82a) 

Ml 

e  s 

to 

V* 

a)  expressed  as  a  fraction  of  Wp- 

(82b) 

S  _  0)’  , 

(o’  expressed  as  a  fraction  of  u>2. 

(82c) 

o)2 


The  generalized  expressions  then  become 

1 

EN^$ss^r^  *  2cI  f  (1_x)  "'"'friqx^  cos  2iIqxr  dx  (83) 

o  4 


C0VNt<!)ss(r)’  *ss(8)]  = 


4q 


2  < 


o;1(l-x)  ^(l-xICr-fs)  cosn  (r.s)  dx 

2llqx  Ilq  (l-x)(r+s) 


,  sin  20 qx  sin  n_q(l-x)  (r-s).  Haxfr+s')  dx 

/  (1_x)  2IIqx  Ilq  ( 1— x)  (r-s)  C0S  qX(r+S;  dX 


(84) 


VarN[$ss(r)]  -  (E^ (r) ] }2 


+  4q‘ 


/  (1-x)  s*n  21Iqx  sin  2IIqr(l-x)  dx 
,  2IIqx  2Hqr(l-x) 


(85) 
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with  a  corresponding  expression  for  the  spectral  variance. 

A  Fortran  IV  computer  program  was  written  to  calculate  and  plot  the  normalized 
results  derived  above.  Because  of  the  symmetry  properties  of  Cov^[$ss(r) ,  $sg(s)], 

it  is  necessary  to  study  the  behavior  of  the  function  only  over  the  upper  right- 
hand  quadrant  of  the  r-s  plane.  Similarly,  since  %r*ss  (r)J  is  an  even 
function  of  r,  it  is  presented  only  for  r  positive.  The  reader  will  recall  that 
this  data  applies  to  the  case  where  an  input  noise  process  with  a  rectangular  low 
pass  spectrum,  as  in  Figure  6a,  is  gated  by  a  rectangular  pulse  of  duration  T,  as 
in  Figure  4.  The  most  important  parameter  of  the  situation  is  q,  defined  to  equal 
the  number  of  periods  of  the  highest  noise  frequency  contained  in  the  pulse 
duration  T  -  the  higher  q  becomes,  the  longer  becomes  the  "piece"  of  the  input 
process  gated  by  the  pulse.  The  results  are  plotted  against  frequency  variables 
expressed  as  a  percentage  of  the  highest  noise  frequency  io2 .  Hopefully,  this 
approach  will  divorce  the  data  from  particular  parameter  values. 

Figure  7  shows  a  plot  of  EM[$  (r)]  for  this  case  with  q  =  10.0.  Because 

JN  S  S 

q  is  appreciably  larger  than  unity,  reciprocal  spreading  arguments  show  that  the 
energy  density  spectrum  of  e(t)  is  considerably  narrower  than  the  power  spectrum 
of  n(t).  Thus,  the  convolution  implied  in  the  calculation  of  yields 

a  spectrum  of  approximately  the  same  shape  and  width  as  that  of  n(t).  (The  width 

of  E[$sg(oi)]  is  therefore  approximately  to2 ,  and  that  of  E^[  $sg  (r)  ]  about 

a) 2  /  =  1  in  units  of  the  variable  r.)  The  normalized  expected  spectrum  is  seen 

/(i)p 

to  be  roughly  rectangular  with  amplitude  approximately  1/2  and  half-width  unity. 
This  checks  the  normalization  which  was  intended  to  set  the  area  under  the  curve 
equal  to  unity.  The  "tail"  of  the  spectrum,  extending  beyond  r  =  1.0,  is  the  most 
obvious  effect  of  the  gating  process  on  the  input  spectrum. 


Figure  8  shows  a  plot  of  the  standard  deviation  of  <t>ss(r)  as  a  fraction  of 


the  mean  of  $  (r) ,  i.e., 

S  6 


i(r)  - 


/arN[»ss(r) 3 

F.N[*ss<r)] 
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As  noted  above,  this  ratio  must  always  be  greater  than  unity,  and  near  the 
origin  achieves  its  maximum  value  of  v2.  This  is  because  it  can  be  shown  from 
equation  (75)  and  the  first  form  in  equation  (78)  with  w  =  0  that 


VarN[4,ss(0)]  =  2{EN[$ss(0)]}2  (89) 

The  fact  that  for  large  values  of  r,  a^(r)  ^  1  implies  that 

VarN($ss(r)]  ^  (EN[$ss(r)]}2  for  r»0  (90) 

or  alternatively,  that  l2£(m,w)  of  equation  (78)  makes  negligible  contribution 
to  the  variance  for  o)>>0. 

A  three-dimensional  rendering  of  the  normalized  spectral  covariance  surface 
for  the  upper  right-hand  quadrant  of  the  r-s  plane  is  shown  in  Figure  9.  It  is 
seen  that  the  surface  has  the  form  of  a  sharp,  narrow  ridge  symmetrical  about  the 
plane  r  =  s,  running  from  the  origin  to  the  point  r  =  s  =  1  and  then  tapering 
off.  From  equation  (77),  we  know  that  the  height  of  the  surface  above  the  line 
r  =  s  is  a  measure  of  Var^[<J>ss(r)] .  Accordingly,  the  height  of  the  surface  at 

the  origin  is  roughly  twice  that  found  elsewhere  along  the  line  r  =  s,  in  agreement 
with  equations  (89)  and  (90) . 

As  an  aid  to  relating  the  height  and  length  of  the  covariance  ridge  to  the 
parameters  of  the  random  transient,  let  us  examine  the  r-s  piano  itself  as  shown 
in  Figure  10.  As  noted  several  times  before,  the  covariance  sutrace  is  symmetrical 
about  the  plane  defined  by  the  relationship  r  =  s.  Consider  any  line  in  the 
r-s  plane  parallel  to  that  representing  r  =  s.  The  equation  cf  such  a  line  is 

s  =  a  +  r 

OR  (91) 


s  -  r  =  a 

where  a  is  the  intercept  on  the  s  axis.  Evidently,  along  such  paths,  s  and  r 
have  a  constant  difference  equal  to  a.  In  the  present  case,  such  lines  represent 
the  loci  of  points  for  which  the  arguments  of  the  covariance  function  have  a  given 
(normalized)  difference  equal  to  a.  Similarly,  lines  lying  n -rpendicula:  to  the 
symmetry  axis  have  equations  of  the  >  orm 

s  +  r  =  b  (92) 

and  on  such  loci,  the  sum  of  r  and  s  is  a  constant.  In  particular,  since 
(s  +  r)/2  is  the  arithmetic  average  of  the  two  normalized  frequency  arguments, 
such  a  line  represents  the  loci  nf  points  for  which  the  center  frequency  of  the 
two  is  constant  and  equal  to  b/2.  Furthermore,  it  can  be  shown  that  a  translation 
of  magnitude  d  parallel  to  the  direction  of  the  line  r  =  s  produces  a  change 
in  the  center  frequency  equal  to  d//2~  while  leaving  the  center  frequency  unchanged. 


NORMALIZED  SPECTRAL  COVARIANCE  SURFACE  FOR  A  RECTANGULAR  PULSE 
OF  LOW-PASS  BROAD  BAND  NOISF  WITH  q  =  10.0. 


FIG.  10  IMPORTANT  LINES  IN  THE  r  -  s  PLANE 


FIG .  1 1  NORMALIZED  SPECTRAL  VARIANCE  FOR  A  RECTANGULAR  PULSE  OF 
LOW-PASS  BROAD  BAND  NOISE  WITH  q  =  10.0. 
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'Returning  now  to  the  case  portrayed  in  Figure  9,  Figure  11  shows  essentially 
a  section  through  the  covariance  surface  along  the  plane  of  symmetry  r  =  s.  This 
can  be  interpreted  as  a  graph  of  VarN[$sg(r)]  plotted  along  an  axis  scaled  by  a 

factor  of  JT.  Along  the  line  r  =  s,  the  covariance  arguments  have  zero  differ¬ 
ence,  and  the  covariance  reduces  immediately  to  the  variance.  From  another  point 
of  view,  one  obtains  a  section  through  the  plane  of  symmetry  by  graphing  the 
variance  as  a  function  of  r,  and  then  "stretching"  the  graph  along  the  r  axis 
by  a  factor  of  J2~.  In  Figure  12  are  shown  two  plots  of  CovN[$ss(r) ,  $ss  (s)] 

where  a  constant  center  frequency  has  been  chosen  as  a  parameter  and  the  covariance 
plotted  as  a  function  of  the  difference  between  the  two  arguments.  In  other  words, 
we  portray  Cov^(r?  s)  where 


r  =  c  +  d/2  (93a) 

and 

s  =  c  -  d/2  (93b) 


This  is  an  even  function  of  d,  so  only  half  of  each  plot  has  been  presented.  From 
the  discussion  immediately  above,  it  becomes  obvious  that  such  graphs  represent 
cross-sections  of  the  covariance  surface  perpendicular  to  the  plane  of  symmetry 
and  intersecting  the  latter  above  the  point  r  =  s  =  c.  The  normalized  frequency 
difference  d  represents  a  translation  perpendicular  to  the  plane  of  symmetry  of 
length  d/^2  .  The  family  of  graphs  thus  represents  a  series  of  cross-sections  of 
the  surface,  taken  perpendicular  to  the  axis  of  the  "ridge". 

In  Figures  13  through  17,  a  similar  series  of  graphs  are  shown  for  the  case 
of  a  rectangular  pulse  of  low-pass  broad  band  noise  with  q  =  1.0.  It  is  seen  that 
the  covariance  surface  now  appears  as  a  relatively  smoothly  sloping  "pile"  which 
decreases  gradually  in  all  directions  with  distance  from  the  origin.  Recalling  the 
definition  of  the  parameter  q,  the  difference  between  the  two  random  transient 
examples  treated  thus  far  is  that  in  the  former  case,  where  q  =  10.0,  a  relatively 
large  number  of  cycles  of  the  highest  frequency  present  in  n(t)  are  found  in  pulse 
length  T,  whereas  for  q  =  1.0,  only  a  single  cycle  is  possible.  The  q  =  10.0  case 
essentially  presents  a  longer  "piece"  of  the  noise  process  n(t). 

Some  general  characteristics  of  the  spectra  of  rectangular  noise  pulses  can  be 
immediately  discerned  by  studying  Figures  7  through  17.  It  turns  out  that  an  inter¬ 
esting  parameter  of  the  random  transient  spectrum  is  its  effective  width,  essentially 
defined  as  that  positive  frequency  spread  for  which  E[$ss(w)]  is  significantly 
different  from  zero.  It  will  be  recalled  that  this  expected  value  can  be  expressed 
as  the  convolution  of  the  energy  density  spectrum  of  the  envelope  function  and  the 
power  density  spectrum  of  the  underlying  noise  process.  Roughly  speaking,  the  width 
of  the  convolution  of  two  even  functions  is  equal  to  the  sum  of  the  widths  of  the 
functions  themselves.  The  justification  of  this  approximation  is  discussed  in 
Appendix  B.  If  one  spectrum  is  significantly  wider  than  the  other,  the  width  of  the 
convolution  of  the  two  will  essentially  be  that  of  the  former.  In  the  present  case, 
the  positive  frequency  width  of  the  (double-sided)  noise  process  spectrum  is  simply 
C02  (or  1  in  terms  of  the  normalized  frequency  variables).  The  energy  density 

spectrum  of  a  rectangular  pulse  of  length  T  and  unit  amplitude  is  given  by 
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q- 


( RELATIVE  MAGNITUDE:  1 .0  HERE 

correspond:  to  actual  value 

OF  .247) 

c  =  (r  +  s)/2,  THE  CENTER  FREQUENCY 
d  =  r  -  s,  THE  DIFFERENCE  FREQUENCY 
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FIG.  12  SECTIONS  THROUGH  THE  SPECTRAL  COVARIANCE  SURFACE  FOR  A 
RECTANGULAR  PULSE  OF  LOW-PASS  BROAD  BAND  NOISE  WITH 
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FIG.  17  SECTIONS  THROUGH  THE  SPECTRAL  COVARIANCE  SURFACE  FOR  A  RECTANGULAR 
PULSE  OF  LOW-PASS  BROAD  BAND  NOISE  WITH  q  =  1 .0  . 


NOLTR  68-124 


4*  (w) 

ee 


li 

2n 


sin2 


o)T 

2 


9 


-00  <  (|)  <  00 


(94) 


as  was  shown  In  reference  (a),  page  27.  The  first  zero  for  to  >  0  occurs  when 
the  argument  of  the  sine  in  the  numerator  equals  H  ,  or  when 


2 n 

T 


(95) 


This  point  may  be  taken  to  define  the  positive  frequency  width  of  the  double-sided 
energy  density  spectrum  of  the  envelope,  and  hence  by  the  arguments  presented 
above,  the  effective  positive  frequency  width  of  the  spectrum  of  a  rectangular 
pulse  of  low-pass  broad  band  noise  can  be  said  to  be 


W  =  u2 


(96) 


or  when  normalized  to  the  same  basis  as  the  frequency  variables  r  and  s, 

u2  211 

%  *  —  +  •  1  + 17 1  <97> 

where  q  is  defined  in  equation  (82a).  When  q  =  10,  as  for  the  case  portrayed 
in  Figure  7,  WN  *  1.1,  whereas  for  q  =  1.0,  as  in  Figure  13,  WN  =  2.0.  The 

size  of  the  effective  spectral  width  is  mirrored  directly  in  the  form  of  the 
spectral  covariance  function.  It  is  seen  empirically  that  the  spectral  covariance, 
centered  above  the  line  r  =«  s,  falls  to  zero  beyond  the  point  r  “  s  =  W^.  In 

other  words,  the  covariance  between  the  values  of  the  spectrum  at  two  points  will 
be  small  if  either  argument  lies  outside  the  region  of  significantly  large  expected 
value.  This  is  well  demonstrated  in  the  q  «*  10.0  case  where  the  covariance  surface 
falls  sharply  away  for  either  r  or  s  >  1  (see  Figures  9  and  11)  and  in  the 
q  ■  1.0  case  where  the  surface  is  essentially  of  zero  height  for  either  r  or  s  >  2 
(see  Figures  15  and  16). 

In  seeking  the  width  of  the  covariance  "mound"  in  a  direction  perpendicular 
to  its  axis,  one  sees  from  Figures  12  and  17  that  it  is  on  the  order  of  2/q.  This 
implies  Immediately  that  the  longer  the  random  transient  becomes  with  respect  to 
the  period  of  the  highest  noise  frequency,  the  narrower  becomes  the  spectral  co- 
variance  function  about  its  plane  of  symmetry.  Thus,  when  q  =  10.0,  the  width 
of  the  ridge  is  considerably  less  (by  a  factor  of  10)  than  the  width  when  q  =  1.0. 
The  peak  at  the  origin,  which  arises  as  a  consequence  of  equation  (89)  dies  off  for 
r,  s  >  1/q.  This  is  especially  clear  in  Figure  9  for  the  q  »  10.0  case. 

In  general,  the  volume  lying  under  the  covariance  surface  can  be  said  to  be 
concentrated  over  the  r  -  s  plane  area  shown  in  Figure  18.  This  region  is 
bounded  by  the  lines : 
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r 

ii 

(98a) 

s 

ii 

(98b) 

s 

=  r  +  1/q 

(98c) 

s 

=  r  -  1/q 

(98d) 

The  peak  at  the  origin  lies  roughly  above 

the  area  bounded  by  the  lines 

r  =  1/q 

and  s  =  1/q. 

When  q>>l?  the  different  sections  of  the  covariance  surface  become  separate 
and  well-defined.  This  is  shown  particularly  well  by  Figure  9.  Near  the  origin, 
one  finds  a  high  narrow  peak  over  an  area  roughly  1/q  on  a  side.  The  largest 
section  of  the  surface  is  the  narrow  ridge  that  runs  diagonally  from  the  peak  at 
the  origin  to  the  vicinity  of  the  point  r  =  s  =  1.0,  where  it  decreases  sharply 
to  insignificance.  The  height  of  the  ridge  is  fairly  constant  along  its  length 
and  equal  to  half  that  of  the  peak  at  the  origin.  This  shows  that  for  large  q,  the 
values  of  the  empirical  spectrum  at  two  arguments  do  not  co-vary  strongly  unless 
the  arguments  are  close  together,  and  conversely  that  spectral  amplitudes  for  widely 
differing  arguments  are  practically  linearly  independent. 

When  q  is  on  the  order  of  unity,  1/q  is  of  the  same  order  as  WN,  and  the 
various  regions  of  the  covariance  surface  fuse  together  such  that  no  separate 
identifications  can  be  made.  In  this  case,  the  surface  becomes  much  more  smoothly 
varying  and  spreads  away  from  the  plane  r  =  s,  as  was  seen  in  Figure  15.  This 
implies  that  significant  covariance  exists  between  widely  separated  points  on  the 
spectrum  and  that  a  high  degree  of  dependency  is  present  between  far  removed  sections 
of  the  spectrum. 

The  areas  of  non-negligible  covariance  delineated  in  Figure  18  agree  well  with 
intuition.  Certainly  one  would  expect  a  strong  covariance  between  the  spectral 
amplitudes  of  two  points  only  when  the  expected  spectral  amplitudes  themselves  have 
appreciable  magnitude  and  can  be  said  to  arise  from  something  more  than  chance.  The 
width  of  the  ridge  essentially  indicates  that  two  spectral  values  do  not  co-vary  if 
their  normalized  arguments  differ  by  more  than  Ar  =  1/q.  There  is  a  well  known 
approximate  frequency  resolution  theorem  which  states  that  for  any  random  time 
function  of  duration  T,  energy  density  spectrum  values  are  statistically  independent 
if  their  arguments  have  a  difference  of  at  least  Au  *  2H/T.  (See  for  example 
reference  (f),  section  6.1.1).  Since  by  definition 
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and  in  terms  of  the  unnormalized  frequency  variables. 


Am 


m2  Ar 


2H 

T 


(100) 


Thusf  the  form  of  the  covariance  surface  bears  out  the  prediction  of  this  particular 
frequency  resolution  theorem. 

B .  Rectangular  Pulse  of  Narrow  Band  Noise 

We  now  turn  to  consider  spectral  covariance  functions  for  narrow  band 
random  transients..  again  assuming  a  rectangular  envelope  function.  Expressions  for 
this  noise  spectrum  and  the  corresponding  autocorrelation  function  have  already 
been  given  in  equations  (71)  and  (72).  Again  the  various  statistical  measures  will 
be  presented  in  normalized  versions  in  which 


m0T 

q  =  — ,  the  number  of  periods  of  the  noise  center  frequency 

contained  in  a  pulse  length.  (101a) 


r  =  -  ,  a)  expressed  as  a  fraction  of  w  . 

w0  0 


(101b) 


S  =  0) 


,  oj*  expressed  as  a  fraction  of  o>  . 


(101c) 


Am 


O) 


the  noise  bandwidth  as  a  fraction  of  the  center  (lOld) 

frequency. 


The  resulting  expressions,  obtained  precisely  as  in  the  broad  band  case,  become: 

1 

EmU  (r)]  =  2q  J  (1-x)  -l)  -^Z'X  cos  2Jlqx  cos  2Bqrx  dx  (102) 

ss  liqzx 


such  that 

00 

/  EN[*sg(r)]dr  -  1 

^00 


and 

®t ♦„<«>]  -  %t«ss(r  ■  ir” 

0 


(103) 


39 


NOLTR  68-124 


Covt*gsCr),  *SSC«0] 


4q2 


/( 1-x) 

Ilqzx 


cos 


on nv  sin  Ilq(l-x)  (r+s) 
q  nq (1-x) (r+s) 


/  (1-x)  sl?.n.az2L 
,  Ilqzx 


cos 


„„  sin  Ilq(l-x)  Cr-s) 
2II(1X  - n'q(l-xt(r-sT 


cos  Ilqx(r-s)  dx 


cos  Hqx(r+s)  dx 


(104) 


with 


Cov[  $  (oj)  ,  <J>  (a)' )  ] 

Sb  S  S 


16N  2T2  Cov[$  (r  =  — ),  4>_(r  =  7-)] 


ss 


co_ 


ss 


(105) 


Two  examples  of  the  rectangular  pulse  of  narrow  band  noise  are  shown  in  the 
figures  to  display  the  characteristics  of  the  corresponding  spectral  covariance 
surfaces.  In  Figure  19,  is  shown  the  spectral  expected  value  for  a  narrow-band 
random  transient  with  q  =  50.0  and  z  =  0.1.  This  is  to  say  that  the  length  of 
the  transient  is  such  that  50  cycles  of  the  center  frequency  can  be  contained 
therein,  and  that  the  bandwidth  of  the  underlying  noise  process  is  1/10  the  center 
frequency.  Using  spectral  convolution  arguments  as  before,  the  positive  frequency 
spectral  width  will  be  on  the  order  of 


W  = 


411 


2Aw 


(106) 


and  normalization  with  w  as  the  basis  yields 

0 


Wv 


_2_ 

q 


+  z 


(107) 


In  the  present  example,  =  .14,  as  is  borne  out  by  the  graph  of  Figure  19. 

The  spectral  covariance  surface  for  this  case  is  rendered  three-dimensionally  in 
Figure  21,  and  in  section  in  Figure  22.  A  plot  of  the  normalized  spectral  vari¬ 
ance  is  shown  in  Figure  20.  It  will  be  recalled  that  such  a  graph  can  be  interpreted 
as  yielding  a  section  of  the  covariance  surface  through  the  plane  of  symmetry. 

Because  of  the  relative  compactness  of  E[<t  (w)],  the  covariance  surface  is  a 

short,  siharply  peaked  ridge  symmetrically  placed  about  the  plane  r  =  s  and 
concentrated  between  r  =  s  =  1  -  z  and  r  =  s  =  1  +  z.  As  before,  the  sections 
of  Figure  22  are  taken  at  right  angles  to  the  axis  of  the  ridge. 

The  corresponding  set  of  figures  is  shown  for  the  case  where  q  =  10.0  and 
z  *  0.01.  Here  the  bandwidth  of  n(t)  is  very  small  compared  with  the  center 


FIG.  19  NORMALIZED  SPECTRAL  EXPECTED  VALUE  FOR  A  RECTANGULAR  BURST 
OF  NARROW  B^N  D  NOISE  WITH  q  =  50.0  AND  z  =  0. 1 . 

q  =  50.0  z  =  0.1 


FIG. 20  NORMALIZED  SPECTRAL  VARIANCE  FOR  A  RECTANGULAR  BURST  OF 
NARROW  BAND  NOISE  WITH  q  =  50.0  AND  z  =  0. 1 . 
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CovN  !X(r)'  $ss(s)] 


(RELATIVE  MAGNITUDE:  1 .0  HERE 
CORRESPONDS  TO  ACTUAL  VALUE 
OF  23.00.) 


c  =  (r  +  s)/2,  THE  CENTER  FREQUENCY 
d  =  r  -  s,  THE  DIFFERENCE  FREQUENCY 
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o 


FIG. 22  SECTIONS  THROUGH  THE  SPECTRAL  COVARIANCE  SURFACE  FOR  A 
RECTANGULAR  PULSE  OF  NARROW  BAND  NOISE  WITH  q  =  50.0 
AND  z  =  0. 1 . 
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frequency,  and  also  small  compared  to  the  spectral  width  of  the  envelope  function. 
As  a  result,  the  spectral  expected  value  (Figure  23)  is  practically  identical  with 
that  of  the  envelope  (which  has  the  familiar  sin2x/x2  form),  and  in  turn  the  co- 
variance  peak  is  more  widely  spread  and  smoothly  rounded  around  the  vicinity 
r  =»  s  =  1  (see  Figures  23,  25,  and  26).  Here  the  spectral  width  is  approximately 
equal  to  2/q  since  z  is  relatively  so  small. 

In  general,  it  is  found  for  the  rectangular  narrow  band  case  that  the 
covariance  peak  is  found  over  an  area  delineated  in  the  r  -  s  plane  by  the  lines 


r  =  1  +  Wn/2 
r  =  1  -  Wm/2 


(108a) 

(108b) 


s  =  1  +  Wn/2 

s  -  1  -  Wn/2 

s  =  r  +  1/q 

s  =  r  -  1/q 


(108c) 
(108d) 
(108e) 
(108f ) 


as  in  Figure  27.  Again,  this  representation  agrees  with  the  interpretation 
offered  for  Figure  18  and  the  rectangular  broad  band  case. 

C.  Random  Transient  Ensembles  with  Decaying  Exponential  Envelopes 


A  number  of  examples  of  EN[<J>ss(r)J  and  CovN[<i>ss(r) ,  <J>sg(s)]  have  been 
computed  and  graphed  for  random  transient  ensembles  generated  with  decaying 
exponential  envelopes.  In  these  studies,  the  parameter  q  has  been  defined, 
analogously  to  equations  (82a)  and  (101a),  with  T  now  taken  to  represent  the 
time  constant  of  the  exponential  decay,  q  here  represents,  then,  the  number  of 
periods  of  a  typical  noise  frequency  contained  in  the  time  constant.  For  a  given 
underlying  noise  process,  the  various  spectral  statistics  for  the  exponential 
envelope  cases  do  not  differ  in  any  gross  manner  from  their  counterparts  in  the 
rectangular  envelope  cases  with  identical  q(  i.e.,  from  those  arising  from 
rectangular  random  transients  whose  envelope  length  is  equal  to  the  exponential 
time  constant).  The  only  appreciable  difference  observed  is  that  the  spectral 
width  is  slightly  narrower  in  cases  where  the  spectral  width  of  the  envelope  is  a 
significant  part  of  the  total,  and  that  the  transverse  width  of  the  spectral  co- 
variance  ridge  is  somewhat  reduced,  although  it  is  still  well  approximated  by  2/q. 
These  effects  are  demonstrated  in  Figures  28  ard  29  which  compare  the  results  found 
for  exponential  and  rectangular  envelope  narrow  band  random  transients  with  q  -  10.0 
and  z  =  0.1.  These  changes  are  not  unexpected,  since  the  energy  density  spectrum 
of  a  decaying  exponential  pulse  is  significantly  narrower  than  that  of  a  rectangular 
pulse  whose  length  is  set  equal  to  the  time  constant.  This  is  shown  in  Figure  30. 


One  may  say,  then,  that  the  approximate  delineation  of  the  region  of 
appreciable  covariance  set  forth  for  the  rectangular  transients  (Figures  18  and  27) 
applies,  for  most  practical  purposes,  to  the  exponential  envelope  cases  also.  It 
may  safely  be  concluded  that  the  detailed  shape  of  the  envelope  function  is 
relatively  unimportant  for  the  form  of  the  spectral  statistics,  and  that  the  envelope 
duration  -  or  some  consistent  measure  thereof  -  is  the  most  important  envelope 
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FIG. 23  NORMALIZED  SPECTRAL  EXPECTED  VALUE  FOR  A  RECTANGULAR 
BURST  OF  NARROW  BAND  NOISE  WITH  q  =  10.0  AND  z  =  0.01 .  . 


FIG.  24  NORMALIZED  SPECTRAL  VARIANCE  FOR  A  RECTANGULAR  BURST 
OF  NARROW  BAND  NOISE  WITH  q  =  10.0  AND  z  =  0.01 . 


NORMALIZED  SPECTRAL  COVARIANCE  SURFACE  FOR  A  RECTANGULAR  PULSE 
OF  NARROW  BAND  NOISE  WITH  q  =  10.0  AND  z  =  0.01 . 
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feature  for  the  determination  of  the  spectral  character  of  the  random  transient 
ensemble..  This  is  particularly  true  when  the  transient  duration  is  such  that  it 
contains  many  cycles  of  a  typical  noise  frequency,  i.e.,  when  the  parameter  q  is 
large. 

This  concludes  our  discussion  of  random  transient  spectral  covariance  surfaces. 
It  is  hoped  that  the  reader  has  been  given  sufficient  grasp  of  the  concepts  involved 
to  permit  him  to  continue  on  his  own  to  explore  special  cases  of  interest. 
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Chapter  IV 

IDEAL  BAND  PASS  FILTERING  OF  RANDOM  TPANSIENT  WAVEFORMS 


In  this  chapter  will  be  studied  the  problem  of  passing  an  ensemble  of  random 
transients  through  an  idealized  narrow  band  filter  whose  system  function  is  given 
by  the  expression 

H(w)  =  1,  cuc - -  <  ja)J  l(iic  +  -|-  (109) 

=  0,  elsewhere 


where  w is  the  center  frequency  of  the  filter,  and  B  is  its  bandwidth,  both 
in  radians/second.  B(w)  is  portrayed  in  Figure  31. 

A  fow  words  are  in  order  about  this  filter  function.  First  we  note  that  such 
a  device  is  physically  non-realizable.  The  impulse  response  h(t)  is  given  by  the 
Fourier  transform  of  H(w) : 


h(t) 


00  , 

J  H(-ju)  eJwt  dw 

—rc 


wc+B/2 

^1 _  /  ,  cos  u)£  do) 

n  «2-b/2 


COS  u)  t» 


—  00 


t  <  00 


(110) 


and  is  non-zero  for  t  0,  implying  that  the  device  responds  before  the  arrival  of 
a  stimulus.  This  is  a  physical  impossibility  (at  least  in  this  world) .  Since  the 
filter  of  ecuation  (10Q)  is  unrealizable,  one  might  wonder  why  it  will  be  studied 
here.  Of  course,  the  honest  answer  is  that  it  is  mathematically  tractable,  but  it 
should  be  pointed  out  that  if  an  arbitrarily  long  time  delay  is  allowed  between  the 
input  and  output  (essentially  by  relaxing  the  specification  to  permit  non-zero 
phase),  such  a  filter  function  can  be  realizably  approximated  to  any  accuracy  (see, 
for  example,  reference  (g),  pp  240-241),  At  the  same  time,  the  filter  function  of 
equation  (109)  is  a  convenient  mathematical  abstraction  which  provides  a  means  for 
describing  the  band-wise  content  of  a  random  transient  spectrum  and  particularly 
the  degree  of  variability  to  be  expected  in  band-pass  filtering.  This  will  provide 
a  strong  indication  of  the  magnitude  of  these  effects  in  real  world  cases  with 
realizable  filter  functions. 
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FIG. 31  IDEAL  NARROW-BAND  FILTER  FUNCTION. 


FIG. 32  VAR  CF(Hfl  PORTRAYED  AS  THE  VOLUME  OF  A  PRISM  LYING  UNDER  THE 
SPECTRAL  COVARIANCE  SURFACE. 
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With  this  justification,  we  shall  proceed.  From  equations  (22)  and  (109),  one 
may  immediately  find  that 


-wc+B/2  -c 


F(H)  = 


to  „+B/2 


f  4>ss(<o)  dto  +  f  <J>ss(to)  dw 
“c-B/2  toc-B/2 


(111) 


and  thus  F(H)  is  indeed  an  expression  for  the  area  of  4>  (to)  lying  between  the 

s  s 


frequencies  (to  -  B/2)  and  (to  +  B/2).  Since  $  (to)  is  an  even  function  of  to 
c  c  ss 


toc+B/2 


F(H)  -  2  /  $  (u)dw 

toc-B/2  Sa 


and  further, 


o>c+B/2 


E[F(H)  ]  =  2  /  E[$ss(to)]dto 

Similarily,  from  equation  (30) 


toc-B/2 


~<oc+B/2  -toc+B/2 


Var[F(H) ]  =  f  f  Cov[<Dss(to),  4>ss(to')]  dm  dto* 


-wc-B/2  -coc-B/2 


-wc+B/2 

»c+B/2 
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f 

-toc-B/2 

mc-B/2 

toc-!-B/2 

-mc+B/2 
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/  i 

toc-B/2 

-toc-B/2 

toc+B/2 

toc+B/2 

/ 

/  l 

to  -B/2 

to  -B/2 

rss'w/»  ss 


f  Cov  [ <J>SS  (to)  ,  <I>ss(a)')]  dto  dto' 


ss  ’  ss 


(112) 


(113) 


(114) 
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By  the  symmetry  properties  of  the  covariance  function,  this  becomes 

a>c+B/2  mc+B/2 

Var[E(H) ]  =  4  /  /  Cov[$ss(io),  4>sg (to* ) ]  da)  dm'  (115) 

o)c~B/2  (joc~B/2 

Thus,  the  expected  value  of  the  output  energy  is  given  by  twice  the  area  under 
the  curve  of  E[$ss(<d)]  between  w  =  wc  -  B/2  and  w  =  u>c  +  B/2  whereas  the 

variance  of  the  output  energy  is  found  to  be  equal  to  four  times  the  volume  of  a 
prism  bounded  above  by  the  covariance  surface  and  lying  over  a  square  in  the 
(o  —  w’  plane  bounded  by  the  lines 


0)  =  0)c  +  B/2 

(116a) 

to1  =  o>c  +  B/2 

(116b) 

as  shown  in  Figure  32.  Having  arrived  at  an  expression  for  Cov[4>ss  (w)  ,  $ss(u)’)] 
using  the  methods  of  the  last  chapter,  it  is  straightforward,  at  least  conceptually, 
to  compute  E[F(H)]  and  Var[F(H)]. 

In  practice,  an  accurate  computation  of  the  output  energy  expected  value  and 
variance  is  an  exceedingly  tedious  business.  Several  attempts  have  been  made  to 
perform  the  calculation  by  digital  computer,  using  multiple  numerical  integration. 

On  the  IBM-7090  at  NOL,  the  computation  time  for  a  single  point  on  a  typical  graph 
of  Var  [F(H)]  was  on  the  order  of  30  minutes!  It  appears,  then,  that  an  accurate 
computer  calculation  of  these  measures  is  completely  out  of  the  question,  and  one 
is  forced  to  the  use  of  approximations  which,  while  lacking  accuracy,  nonetheless 
indicate  and  preserve  the  magnitudes  of  the  important  trends  in  the  various 
statistical  functions.  An  example  of  such  an  approximate  method  will  be  detailed 
now. 


For  our  example,  we  shall  examine  the  ideal  frequency  filtering  of  a  rectangu¬ 
lar  pulse  of  narrow  band  noise,  characterized  by  the  parameters  q  and  z  defined 
previously.  The  spectral  covariance  surface  of  such  a  signal  has  been  treated  in 
the  preceding  chapter,  and  its  features  portrayed  in  Figures  19  through  27.  Of  all 
the  examples  treated  in  this  study,  this  is  the  one  of  most  widespread  interest 
and  generality,  particularly  since  the  results  retain  sufficient  accuracy  to  pre¬ 
dict  the  effects  of  filtering  any  shaped  pulse  of  narrow  band  noise  of  approximately 
the  same  duration  as  the  rectangular  pulse,  in  accordance  with  the  discussion  at 
the  close  of  Chapter  III. 

In  the  discussion  that  follows,  the  normalized  and  parameterized  versions  of 
the  expected  values  and  covariances  will  be  used.  To  start  the  present  approximate 
approach,  the  spectral  expected  value  (of  Figure  19  or  23)  will  be  represented  as 
a  rectangular  function  of  width  centered  at  r  =  1.0.  The  height  of  this 

function  will  be  taken  to  be  the  maximum  value  of  the  expected  value  expression 
computed  from  equation  (102).  This  maximum  will  be  denoted  Mj?,  and  thus 
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EAfWr)l  =  me»  1  -  %/2  1  |r|  £  1  +  Wn/2 


(117) 


=  0,  elsewhere 

where  EA[$ss(r)]  represents  the  approximate  normalized  spectral  expected  value. 


E[F(H)]  is  found  from  the  expression  of  equation  (113): 

ojc+B/2 

E[F(H)]  -  2  /  E[4>ss(w)3  du 

mc-B/2 

Using  also  equation  (103)  and  defining  ioc*  =  ojc/o)0  and  B'  =  B/“  , 

o)c+B/2 

E[F(H)  ]  =  8N0T  /  EN[$ss(r  =  to/co0)  ]  dw 

mc-B/2 

a)c'+B'/2  (118) 

=  8NQa)  T  /  E[*  (r)]  dr 

coc'-B72  SS 

Now  denoting  the  approximate  value  of  E[F(H) ]  by  EA[F (H) ] ,  we  have  that 

«c'+B72 

E[F(H)]  «  Ea[F(H) ]  =  k  /  EA[<J>ss(r)]  dr  (119) 

wc'-B72 

where 

k  =  8N0o)0T  (120) 


The  integral  of  equation  (119)  can  be  visualized  as  proportional  to  the  area  of 
overlap  of  two  rectangular  functions:  one  fixed  -  representing  the  random  transient 
spectrum  (equation  (117)),  and  the  other  sliding  -  representing  the  tunable  filter 
as  the  latter's  center  frequency  is  swept  across  the  frequency  range  of  interest. 
This  is  shown  in  Figure  33.  The  resulting  output  energy  expected  value  express¬ 
ions  can  be  considered  a  function  of  wc',  or  by  defining 


wc  (TUNABLE) 

Ih(,)I2ea  [#n«] 

E.  [f(H)1  ~  AREA  OF  OVERLAP 

FIG. 33  AN  INTERPRETATION  OF  THE  CALCULATION  OF  E*  [F(H)J  AS  THE  AREA  OF 
OVERLAP  OF  TWO  RECTANGULAR  FUNCTIONS. 
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as  the  normalized  frequency  difference  between  the  center  of  the  filter  function 
and  that  of  the  underlying  noise  spectrum,  as  a  function  of  A.  Clearly,  there 
are  two  cases  to  consider: 

A.  For  B'  <  WK, 

W  B  ^ 

Ea[F(H)  ]  =  MEkB\  0  £  [A  |  _<  ^~2 - 


B’  +  %  * 

=  MEk[- - 2 - N  : 


WN  -  B' 


(122) 


=  o,  |a]  >_ 


B’  +  Wv 


B.  For  B’  1  Wjj, 


Ea[F(H)J  = 


MEkWN ’ 


0  <  )A|  < 


_  M  i  r  +  %  i  .  B*  -  WN  .  . 

=  MEk[ - 5 ^ —  |a|]  , - 5-2-1  |a|  1 


B’  +  Wv 


(123) 


=  o, 


B'  +  WN 
2 


These  functions  are  shown  in  Figure  34  and  represent  approximate  expressions  for 
the  expected  value  of  the  filter  output  energy  as  a  function  of  the  filter  center 
frequency  wc?  when  the  input  is  a  rectangular  pulse  of  narrow  band  noise  centered 
at  . 


Corresponding  approximate  expressions  involving  the  spectral  covariance 
function  are  somewhat  harder  to  formulate  and  handle.  For  our  present  needs,  the 
spectral  covariance  surface  can  be  reasonably  well  approximated  by  a  right  hexagonal 
prism  erected  over  a  base  in  the  r  -  s  plane  delineated  by  the  equations  (108)  and 
shown  in  Figure  27.  The  height  of  this  prism  is  taken  to  be  the  maximum  value  of 


the  expression  of  equation  (104),  denoted  Mc< 
the  approximated  surface  is  shown  as  Figure  35. 


A  three-dimensional  rendering  of 


The  variance  of  the  filter  output  energy  is  expressed  by  means  of  equations 
(115)  and  (105): 


V 


B‘-wn 


W  +B1 


CJ  -U) 
c  < 
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t»)c+B/2  o)c+B/2 

Var[F(H)]  =  4  /  /  Covt4>„«, (u) ,  iccCu1)]  dm  dto* 

»c-B/2  »c-B/2 


ujc+B/2  0)c+B/2 

=  64  N  2T2  /  /  CovN[*ss(r  =  jp),  $ss(s  =  p-)]duj  doi* 

ioc-B/2  u)c-B/2  0  ° 


co,c+B,/2  a) * C+B ' / 2 

=  64  N  2T2a)  2  /  / 

w'  -B72  to'  -B’/2 

C  c 


CovN^ss(r)>  $ss(s)]dr  ds 


(124) 


If  we  denote  the  approximate  covariance  surface  as  C°vA^ss^r^’  ^ss^s^’ 

to'  +B'/2  o)’+B,/2 
c  c 

Var[F(H)]^  VarA[F(H)]  =  k2  /  /  Cov.[$  (r) ,  $ss(s)]dr  ds  (125) 

(o'c-B'/2  co ' C“B ' / 2 


where  k  is  again  given  by  equation  (12) .  The  integral  of  equation  (125)  can  be 
interpreted  as  yielding  the  volume  under  the  approximate  covariance  surface  lying 
above  a  square  in  the  r  -  s  plane  bounded  by  the  lines: 


r  =  (otc  +  B'/2 
s  =  co’c  +  B'/2 


(126) 


Note  that  this  "filter  square"  has  the  same  axis  of  symmetry  (r  =  s)  as  the  spectral 
covariance  surface  itself.  A  geometrical  portrayal  of  the  generation  of  the  curve 
of  Var^[F(H)]  versus  (o'c  is  sh'wn  in  Figure  36.  There  we  see  that  as  coc 

changes,  the  filter  square  slides  up  and  down  the  line  r  =  s,  cutting  out  a 
changing  volume  under  the  covariance  surface  proportional  to  the  output  energy 
variance.  In  the  present  approximation,  where  the  surface  is  of  constant  height 
over  a  sharply  delineated  area,  outside  of  which  it  is  identically  zero,  the 
integral  of  equation  (125)  becomes  proportional  to  simply  the  area  of  overlap 
of  the  filter  square  and  the  hexagon  representing  the  base  of  the  covariance 
surface.  The  constant  of  proportionality  between  the  overlap  area  and  the  approxi¬ 
mate  variance  is  the  height  of  the  surface,  taken  here  to  be  Mc,  the  maximum 
normalized  covariance.  To  compute  Var^[F(H)]  as  a  function  of  the  filter  center 
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frequency,  one  lays  out  the  covariance  hexagon  as  a  function  of  the  random  transient 
parameters  (q  and  z)  and  calculates  the  area  of  overlap  as  the  filter  square  glides 
up  and  down  the  45°  line  in  the  r  -  s  plane.  Since  Wj$  =  2  +  z  >  JL  always,  the 

q  q 

geometry  of  the  situation  leads  to  three  distinct  cases  depending  upon  the  normal¬ 
ized  filter  band-width  B': 

A.  0  £  B'  £  1/q: 


Var  [F(H)3  =  M  k2B'2  ,  0  | A  |  <_ — 


WM-B' 


Mck^ 


A2-(B'+Wn)[a|-<-  Hril  <|A|<  %+B.l  (127) 


=  0,  A  > 


b'+wn 


B.  1/q  <  B'  <  WN: 

VarA[F(H)]  =  Mck2 


2B 1  _  1 

q  q7 


,  0  <  A  < 


WN-B> 


Mck2 


B'+Wn-2|a|  1 

-  -  — J 

q  q 


wn-b’ 


<|A|<^-_1_  (128) 


Mck2 


A2-(B'+Wn)  1 A  1  +  -(B  ^-2 


WN+B'  1  <  |  A  |  <  WN'+B' 
2  q  2 


=  0,  I A | 


> 


b'+wn 

2 
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C.  B’  >_  WN: 


VarA[F(H)]  = 


1 

i7  ’ 


0  <  |A|  < 


B'-Wn 

~T~ 


Mck2 


B’+Wn-2|A 


1  b*-wn  B*+Wn  i 

*  2  ilAll  2  q 


(129) 


=  M„k2 


(B*+Wn) I  a 


(B'+WN)2 
+  4 


B'+Wn  !  < |a|<  Bt+WN 
2  "  q  2 


-  o, 


A  |  21 


b’+wn 

2 


The  three  resulting  variance  functions  are  shown  in  Figure  (37). 

In  the  light  of  practical  measurement  and  signal  processing  problems,  a  most 
interesting  quantity  in  the  present  context  is  the  standard  deviation  of  F(H) 
normalized  by  the  mean.  This  statistic  essentially  indicates  the  spread  of  the 
probability  distribution  of  F(H),  and  hence  the  amount  of  variability  to  be 
expected  in  the  filter  output  energy  from  sample  to  sample  of  the  random  transient 
ensemble.  The  approximate  normalized  standard  deviation,  here  considered  a  function 
of  B'  and  A,  is  denoted  R(B',  A)  and  is  given  by 


R(B',  A)  = 


✓VarA[F(H)] 

EA[F(H)] 


(130) 


At  the  moment,  we  are  almost  in  a  position  to  compute  this  function.  All  of  the 
approximate  expressions  for  the  expected  value  derived  above  can  be  written  in 
the  form 


EA[F(H)]  =  MEkf,(A,  B',  q,  z) 


(131) 


Similarly  for  the  approximate  energy  variance, 


VarA[F(H)]  =  Mck2f2(A,  B\  q,  z) 


(132) 
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FIG. 37  APPROXIMATIONS  FOR  THE  VARIANCE  OF  THE  OUTPUT  ENERGY  OF  A 
NARROW  BAND  FILTER  WHEN  EXCITED  BY  A  RECTANGULAR  PULSE  OF 
NARROW  BAND  NOISE. 
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Thus, 


R(B',  a) 


/f  2  (A,  B* ,  q,  z) 
Me  f 1 (A,  B’,  q,  z) 


(133) 


In  the  present  example  (i.e.,  that  of  passing  a  rectangular  pulse  of  narrow 
band  noise  through  the  filter),  it  is  found  empirically  (by  direct  calculations) 
that  the  maximum  value  of  E^($ss(r)3  comes  when  r  =  1.  Similarly,  the  maximum 
value  of  CovN[$ss(r),  $ss(s)]  occurs  when  r  =  s  =  1.  Thus  we  may  write,  from 
equations  (102)  and  (104),  that 


sinfiqzx  0 

Mg  =  2q  /  (l-x)-jj— — -  cos^  2llqxdx 


(134) 


M, 


4q2 


1 


sinllqzx  sin  2IIq(l-x) 

f  (T-x)  nqzx  cos  2llqx  Jhq  (l-x5 


(135) 


0 


> 


It  is  also  found  empirically  that  in  equation  (135),  the  first  squared  integral  is 

negligible  compared  to  the  second.  Thus  MC^M|,  and  >4lc/ME  drops  out  of  the 

ratio  of  equation  (133).  R(B',  A)  may  now  be  computed  from  equations  (122),  (123), 

and  (127)  -  (129).  Again,  there  are  three  cases  to  consider,  depending  upon  the 
relative  magnitude  of  B': 

A.  For  0  <  B'  <  1/q: 


R(B\  A)  =  1, 


0  1  N  1 


B'+Wn 

- 


(136) 


undefined  elsewhere. 

B.  For  1/q  £  B'  <_ 


R(B ' ,  A) 


/2B'q-l 

nR7- 


0  <_  I A  |  £ 


W»-B’ 

•> 
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2>/q(B'-fWN  -  2|A|)  -  1  Wn-B'  <  _  1 

q(B’+WN  -  2 j A 1 )  *  2  —II—  2  q 


(137) 


=  h 


Wn+B' 

2 


Wn+B’ 

2 


undefined  elsewhere 


C.  For  B’  •>  WN: 

R(B’,  A) 


<1WN 


,  o  _<  )a|  £ 


b’-wn 

2 


(138) 


2^(6*^  ~  2 1 A  j )  -  1  ^  B*-Wn  B*+Wn  i 

q(B'+WN  -  2|A|)  2  -  |A|  q 


=  1 


B’+%  1  b’+wn 

"T— N  r2- 


undefined  elsewhere 

The  standard  deviation  ratio  is  shown  in  Figure  38  for  the  three  cases  defined 
above . 


Another  way  of  presenting  the  aame  data  is  to  choose  a  given  filter  center 
frequency  (which  plays  the  role  of  a  parameter)  and  then  to  study  the  output  energy 
as  a  function  of  the  filter  bandwidth  for  a  random  transient  ensemble  of  given 
parameters.  In  the  present  example,  if  we  center  the  filter  on  ie0(i.e.,  wc  =  uQ) , 

A  becomes  zero,  and  the  normalized  standard  deviation  as  a  function  of  B*  becomes 


R(B',  0)  =  1,  0  <_  B'  <  1/q 


/2iF  q-1 

B'q 


1 

q 


<  B' 


1% 


/2Wn  q-1 

WNq  ’ 


B'  >  W 


N 


(139) 


This  function  is  plotted  as  Figure  39. 
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FIG. 38  APPROXIMATIONS  FOR  THE  NORMALIZED  STANDARD  DEVIATION  OF 

THE  ENERGY  OUT  OF  A  NARROW  BAND  FILTER  WHEN  EXCITED  BY  A  REC 
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The  expressions  presented  above  fof  E^[F(H)],  Var^fF(H)],  and  R(B’,  A)  only 

roughly  represent  the  actual  functions,  and  one  naturally  inquires  about  the 
accuracy  of  these  approximations.  It  appears  that  E^[F(H)]  and  Var.[F(H)]  are 
useful  principally  for  their  portrayal  of  the  rough  functional  dependence  of  the 
output  energy  mean  and  variance  on  the  filter  band-width  and  center  frequency 
and  the  random  transient  parameters.  Hence,  a  fairly  good  idea  of  the  effects  of 
parameter  changes  can  be  gained  from  the  treatment  given  above,  but  not  too  much 
reliance  should  be  placed  on  the  exact  functional  form  derived  nor  on  the  absolute 
magnitudes  arrived  at  for  the  energy  mean  and  variance.  For  R(B',  A),  on  the 
other  hand,  which  is  probably  the  function  of  most  interest,  several  of  the  un¬ 
certain  factors  cancel  each  other  out,  and  the  approximation  well  represents  the 
actual  state  of  affairs  in  both  magnitude  and  rough  functional  form.  The  author 
feels  that  the  approximation  presented  for  R(B*,  A)  is  correct  to  within  20%  of 
its  actual  value,  with  higher  accuracies  for  small  A  and  cases  where  q  is  much 
larger  than  unity.  The  basis  for  this  estimate  is  a  comparison  of  the  approximate 
results  with  the  values  found  for  several  points  by  an  accurate  numerical  integra¬ 
tion  (at  an  enormous  expense  in  computer  time).  These  agreed  to  within  10%. 
Naturally,  by  devising  more  elaborate  approximations  for  the  spectral  covariance 
surface,  one  could  achieve  any  degree  of  accuracy  required  in  the  final  result. 

For  a  wide  variety  of  purposes,  the  present  degree  of  complexity  is  adequate  to 
provide  considerable  insight  into  the  frequency  filtering  of  random  transient 
waveforms.  At  the  very  least,  these  results  provide  order  of  magnitude  estimates 
for  the  statistical  uncertainty  of  the  filter  output  energies,  and  they  thus  pave 
the  way  for  the  preparation  of  adequate  measurement  programs  and  monopulse  process¬ 
ing  studies.  If  greater  precision  is  required,  the  methodology  for  achieving  it  is 
clear  from  the  above  treatment. 

What,  then,  have  we  found?  Given  a  random  transient  generated  by  multiplying 
a  narrow  band  Gaussian  process  by  a  square-integrable  deterministic  envelope,  the 
results  gained  here  enable  one  to  estimate  with  reasonable  accuracy  the  amount  of 
variation  from  sample  to  sample  in  the  output  energy  of  narrow  band  filters  to 
whose  inputs  such  random  transients  are  applied.  The  results  presented  have  been 
derived  in  particular  for  rectangular  envelopes,  but  as  the  discussion  at  the  end 
of  the  last  chapter  indicates,  they  apply  also  to  any  envelope  function  of 
approximately  the  same  length,  the  detailed  structure  of  the  envelope  playing 
little  role  in  the  final  outcome. 

The  collected  results  can  be  placed  in  best  perspective  by  considering  the 
inter-relationships  of  the  three  band-widths  which  arise  in  treating  the  problem. 
These  are: 

1.  The  band-width  of  the  underlying  noise  process,  2Aw.  In  its  normalized 
version,  this  has  been  denoted  z  (2  =  ^Au)  ^ 

o 

2.  The  band-width  associated  with  the  envelope  spectrum.  If  the  nominal 
length  of  the  envelope  function  is  T,  this  is  on  the  order  of  4I1/T,  or  when 
normalized,  2/q (since  q  =  w0T) 

W  ‘ 

3.  The  filter  band-width  B,  or  in  its  normalized  version  B'  =  B /wc. 
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As  we  have  noted  above,  the  normalized  approximate  average  spectral  width 
for  a  random  transient  ensemble  is  given  by 

W„  -  2  +  2, 

q 

i.e.,  the  sum  of  the  normalized  noise  and  envelope  spectral  width.  The  r  -  s 
plane  area  of  significant  covariance  changes  its  shape  in  accordance  with  the 
relative  magnitudes  of  2/q  and  z.  If  z  is  small  compared  to  2/q,  Wfj^2/q, 
and  the  transverse  width  of  the  covariance  peak  becomes  comparable  to  its  length, 
as  in  Figure  25  (where  z  =  0.1  and  q  =  10.0).  For  larger  z,  the  width  of  the 
ridge  becomes  narrow  compared  to  its  length,  implying  that  i>ss(o>)  and  ^ggCm*) 
covary  strongly  only  when  oi^co'.  Such  a  case  is  shown  in  Figure  21  (where 
z  =  0.1  and  q  =  50.0).  As  noted  above,  the  variance  of  the  filter  output  energy 
can  be  visualized  as  the  area  of  overlap  between  the  "filter  square",  whose  side 
is  length  B',  and  the  covariance  area  referred  to  above  (see  Figure  36).  In 
situations  where  2/q  <<  z  and  Wjjj^z,  it  is  evident  that  this  variance  de¬ 
creases  as  q  increases,  i.e.,  as  the  transient  lengthens  in  terms  of  a  typical 
noise  period,  since  the  width  of  the  ridge  goes  as  1/q  while  its  length,  propor¬ 
tional  to  Wfl,  remains  constant. 

The  square  of  the  output  energy  mean,  E^2[F(H)],  can  be  interpreted  in  turn 
as  being  proportional  to  the  overlap  area  of  the  filter  square  and  a  square  whose 
side  is  of  length  Wj^,  as  shown  in  Figure  40.  (This  emerges  from  equation  (119) 
and  Figure  33).  Thus,  R2(B',  4)  can  be  interpreted  as  the  ratio  of  two  areas 
overlapped  by  the  filter  square  -  that  overlapped  on  the  covariance  hexagon,  and 
that  overlapped  on  the  "expected  value  square".  (See  also  Figure  40.)  From  this 
viewpoint,  it  becomes  immediately  obvious  what  degree  of  statistical  reliability 
can  be  expected  from  measurements  on  random  transients  represented  in  the  r  -  s 
plane  as  covariance  hexagons.  We  see,  for  example,  that  R(B',  A)  can  never  exceed 
unity  since  the  covariance  hexagon  is  contained  entirely  within  the  expected  value 
square.  Furthermore,  when  B*  _<  1/q.  the  filter  square  can  be  contained  entirely 
within  the  covariance  area,  and  R(B',  A)  must  always  be  unity,  as  shown  by 
equation  (136)  or  Figure  39.  This  immediately  implies  a  very  significant  fact: 
that  it  is  completely  impossible  to  obtain  consistent,  i.e.,  repeatable,  filter 
output  energies  for  an  arbitrarily  narrow  filter  band-width.  For  high  spectral 
resolution,  naturally  B'  should  be  as  narrow  as  possible,  but  this  very  narrowness 
destroys  the  stability  of  the  measurement,  and  a  trade-off  becomes  decidedly  neces¬ 
sary.  An  uncertainty  principle  is  at  work  here:  the  more  narrowly  we  try  to 
isolate  the  amount  of  spectral  energy  in  a  given  frequency  band,  the  less  precise 
the  estimate  becomes  in  a  statistical  sense. 

At  the  cost  of  spectral  resolution,  more  stable  estimates  can  be  produced  by 
increasing  B',  i.e.,  the  band-width  of  the  analyzing  filter.  This  is  precisely 
demonstrated  in  Figure  39  for  the  special  case  of  u)'c  =  1  (A  =  0) .  For  given  q 
and  Wj,j,  however,  there  is  a  lower  limit  to  the  normalized  standard  deviation, 
given  by 


WB'-  = 


/2qWj]-l 


/2qz+3  /n  /2TAo)+3n 

nrr+0  TAm4-?TT 


(140) 
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reached  when  Bl  =  W^.  In  Figure  40,  this  is  when  the  filter  square  becomes  the 
same  size  as  the  expected  value  square.  Note  that  at  this  point  virtually  all 
spectral  resolution  has  been  sacrificed,  and  still  the  spectral  estimate  cannot 
be  made  arbitrarily  precise.  This  is  a  direct  consequence  of  dealing  with  random 
functions  of~  finite  length. 

For  given  q  and  with  B*  <  Wjj  (in  an  attempt  to  salvage  some  spectral  resolu¬ 
tion)  the  minimum  value  of  normalized  standard  deviation  achievable  is  given  by 
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/2Btq-l 

qB* 
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Note  that  from  our  parameter  definitions, 
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Thus,  for  large  B'q  products, 
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which  is  essentially  the  result  found  in  reference  (h) ,  section  5.4.5,  for  the 
normalized  standard  deviation  of  a  spectral  measurement  performed  with  the  band¬ 
width  of  the  analyzing  filter  (B)  small  compared  with  the  band-width  of  the  process, 
when  the  filter  output  is  observed  and  averaged  for  T  seconds.  An  alternative 
approach  for  this  special  case  is  presented  in  Appendix  C,  Equation  (143)  is  a 
theoretical  confirmation  of  the  rule-of-thumb  frequently  quoted  in  spectral  analysis 
that  a  large  band-width-time  product  is  necessary  for  statistical  accuracy. 


The  uncertainty  principle  referred  to  above  may  be  stated  explicitly  as  follows 
if  the  bandwidth  of  the  analyzing  filter  is  less  than  the  effective  spectral  width 
of  the  ensemble,  i.e.,  if  B'  <  WN,  the  minimum  normalized  output  energy  standard 
deviation  is  given  by  equation  (141).  We  have  also  seen  that  for  large  B'q  products 
that 
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which  implies  immediately  that 
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If  T  is  considered  a  constant  of  the  situation,  then  we  may  say  that 
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BR2(B’,  A)  >  ■—  ,  B'  <  WN  and  B'q  >  >  1 


(145) 


that  is,  that  the  product  of  the  analyzing  bandwidth  and  the  square  of  the  normal¬ 
ized  energy  standard  deviation  (the  normalized  energy  variance)  will  always  be 
greater  than  4II/T  (which  may  be ^identified,  as  in  equation  (106),  as  the  effective 
double-sided  spectral  width  of  a  square  pulse  of  length  T) .  R2(B’,  A)  can  be 

interpreted  as  a  measure  of  the  precision  with  which  the  total  energy  in  the  band 
of  interest  can  be  specified,  since  it  represents  the  variance  of  the  probability 
distribution  for  the  energy  in  the  band.  The  smaller  this  variance  becomes,  the 
less  variation  is  found  in  the  band-wise  energy  content  of  different  samples,  and 
thus  the  more  precisely  this  energy  can  be  predicted  before  the  event.  Similarly, 

B  reflects  the  extent  to  which  the  fine  structure  of  the  ensemble  spectrum  can  be 
resolved  by  the  filter,  since  this  is  the  bandwidth  over  which  the  spectrum  will  be 
smoothed  or  averaged.  Naturally,  spectral  details  of  width  B  or  smaller  will  not 
be  distinguished  by  a  filter  of  that  band-width,  but  the  smaller  that  B  becomes, 
the  higher  becomes  the  spectral  resolution  obtainable.  We  of  course  would  like  B 
and  R2(B',  A)  to  be  simultaneously  as  small  as  possible  so  that  spectral  resolu¬ 
tion  and  amplitude  stability  are  both  high.  Equation  (145)  demonstrates  immediately, 
however,  that  it  is  impossible  to  attain  arbitrarily  high  precision  in  both  resolu¬ 
tion  and  amplitude  at  the  same  time  since  the  product  of  the  frequency  and 
amplitude  "uncertainties"  must  always  exceed  411  /T .  Thus,  a  decrease  in  the 
uncertainty  of  the  spectral  amplitude  comes  at  the  expense  of  a  concomitant 
increase  in  the  averaging  band  and  thus  a  decrease  in  spectral  resolution.  The 
single  parameter  that  determines  the  lower  limit  of  the  uncertainty  product  is  the 
effective  length  of  the  envelope  function  T,  and  in  general,  the  simultaneous 
uncertainty  varies  in  inverse  proportion  to  this  quantity. 

It  is  also  seen  above  that  there  is  a  definite  lower  limit  to  the  normalized 
standard  deviation  for  an  ensemble  of  given  parameters,  and  this  is  reached  when 
B*  =  W^.  This  lower  limit  is  precisely  the  standard  deviation  of  tne  total  energy 
of  an  ensemble  member  and,  our  uncertainty  principle  notwithstanding,  cannot  be 
reduced  by  further  sacrifices  of  frequency  resolution,  since  indeed  there  is  none 
left  to  sacrifice.  The  situation  here  is  such  that  the  uncertainty  relation  stated 
is  not  unrestricted,  and  that  not  only  are  we  prevented  from  attaining  arbitrarily 
high  amplitude  precision  and  frequency  resolution  simultaneously,  but  there  is  even 
a  limit  on  the  extent  to  which  the  latter  can  be  traded  for  the  former.  This  limit 
is  given  as  a  function  of  the  ensemble  parameters  in  equation  (140). 

The  foregoing  should  provide  the  interested  reader  with  at  least  some  insight 
into  the  spectral  distribution  of  energy  in  narrow  band  random  transients.  The 
extension  to  the  broad  band  case  is  straightforward  and  leads  to  the  same  basic 
principle  -  that  a  trade-off  between  spectral  resolution  and  statistical  reli¬ 
ability  is  always  necessary,  and  that  there  is  a  certain  lower  limit  to  the 
normalized  energy  standard  deviation  which  is  determined  by  the  parameters  of  the 
random  transient  ensemble,  and  which  cannot  be  avoided  no  matter  what  kind  of 
filtering  is  employed.  In  all  cases,  as  B'  ->  0,  the  filter  output  energy 
approximates  the  so-called  periodogram  for  4>ss  (co) ,  for  which  it  can  be  shown  that 
R(B',  A)  21  1>  making  it  an  inconsistent  estimate  of  |J>sc(m)  in  the  statistical 
sense  (see  reference  (c),  pp  107-108). 
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Chapter  V 

CONCLUSIONS  AND  SUGGESTIONS  FOR  FUTURE  RESEARCH 


The  theory  set  forth  in  the  present  report  is  a  logical  extension  of  the 
methods  outlined  in  its  predecessor  (reference  (a)),  and  the  treatment  of  random 
transient  frequency  filtering  given  here  has  wide  application  in  areas  of  present 
interest.  These  fall  into  two  broad  classes:  spectral  analysis  of  transient 
random  waveforms;  and  monopulse  signal  processing  applications  where  the  effects 
of  sample- to-sample  variations  are  important.  These  will  be  treated  separately. 

One  of  the  central  problems  of  traditional  spectral  analysis  is  the  empirical 
estimation  of  the  power  spectrum  of  a  continuing  stationary  process,  given  a  sample 
record  of  some  finite  length  T.  The  "intuitive"  approach  of  computing  the  squared 
magnitude  of  the  empirical  Fourier  transform  of  the  sample  record  and  dividing  by 
T  provides  the  so-called  periodogram  of  the  record.  Unfortunately,  the  periodo- 
gram  is  a  random  variable  for  every  value  of  w  which  has  at  all  points  a  standard 
deviation  at  least  as  large  as  the  mean,  no  matter  how  long  the  sample  record  is 
taken.  Contrary  to  intuition,  then,  the  periodogram  is  not  a  consistent  estimator 
of  the  power  spectrum  as  T  increases  without  limit.  For  this  reason,  a  number' 
of  alternative  methods  have  been  devised  which  provide  power  spectral  estimates 
averaged  over  a  band  of  frequencies  as  opposed  to  attempting  a  point  estimate. 

In  the  well-known  work  of  Blackman  and  Tukey  (reference  (i)),  these  estimates  are 
obtained  by  weighting  the  empirical  autocorrelation  function,  before  it  is  Fourier 
transformed,  to  yield  band-integrated  power  spectrum  estimates  of  high  statistical 
reliability.  Similarly,  Bendat  and  Piersol  (reference  (h))  provide  a  rough 
description  of  the  use  of  narrow  band  filters  for  empirical  spectrum  measurement. 

In  both  approaches,  it  is  found  that  the  price  paid  for  higher  statistical  reli¬ 
ability  is  longer  sample  lengths,  or  decreased  spectral  resolution,  or  both.  For 
resolving  the  spectral  fine  structure  of  the  process  of  interest,  the  analyzing 
band-width  should  be  as  narrow  as  possible.  Unfortunately,  the  narrower  this 
bandwidth,  the  less  precise  the  estimate  becomes.  For  ease  of  computation  and  data 
handling,  the  sample  length  should  be  as  short  as  possible.  Unfortunately,  the 
shorter  the  sample  taken,  the  larger  becomes  the  standard  deviation  of  the  estimate. 
Indeed,  as  noted  above  at  the  close  of  Chapter  IV,  this  standard  deviation,  as  a 
percentage  of  the  mean,  is  on  the  order  of  I/^bT  ,  and  for  arbitrary  accuracy,  a 
trade-off  between  the  two  factors  is  necessary. 

Many  of  these  same  considerations  apply  to  the  analysis  of  energy  spectra 
for  random  transient  waveforms,  and  indeed  the  estimation  of  a  process  spectrum 
from  a  sample  function  of  length  T  can  be  considered  a  special  case  of  the 
transient  spectral  problem:  that  in  which  a  rectangular  envelope  function  of 
length  T  is  selected,  and  one  seeks  from  the  separate  realizations  of  the 
resulting  random  transient  ensemble  to  estimate  '^(to)  for  the  underlying  noise 
process.  There  is,  however,  a  difference  in  emphasis.  In  the  case  of  random 
transient  spectral  analysis,  one  has  no  control  over  T,  or  the  effective  length 
of  the  signal,  and  is  faced  with  the  task  of  making  the  best  of  the  situation  in 
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deriving  a  meaningful  spectral  description.  In  a  theoretical  approach,  where  there 
are  reasonable  grounds  for  assuming  the  form  of  e(t)  and  ^nn(w),  the  present  theory 
provides  expressions  for  the  expected  value  of  the  random  transient  spectrum  at  a 
point  and  also  its  variance.  As  we  have  seen,  this  variance  is  rather  large,  no 
matter  what  the  effective  length  of  the  transient.  We  have  thus  learned  very  little 
about  the  spectrum  of  a  single  sample  since  the  probability  distribution  for  the 
spectrum  at  a  point  is  so  widely  spread.  One  way  out  is  to  formulate  predictions 
of  the  spectral  energy  lying  in  non-zero  bands,  since  for  these  the  normalized 
standard  deviation  will  be  on  the  order  of  1/BT  ,  where  B  is  the  analyzing  band¬ 
width  (here  in  hertz).  Naturally,  spectral  resolution  is  sacrificed  in  this 
approach,  but  this  is  the  price  that  must  be  paid  for  statistical  consistency.  In 
this  world,  one  never  gets  something  for  nothing. 

In  an  experimental  approach,  where  one  is  given  a  collection  of  realizations 
of  the  random  transient  process,  and  the  problem  at  hand  is  to  determine  the 
spectral  character  with  no  a  priori  knowledge  about  the  source,  filtering  to  an 
arbitrarily  narrow  frequency  band  destroys  the  repeatability  of  the  results  from 
sample  to  sample,  and  an  average  value  may  or  may  not  constitute  a  reasonable  pre¬ 
diction  for  the  result  to  be  found  in  dealing  with  an  individual  sample.  Probably 
the  best  course  of  action  here  is  to  calculate  the  empirical  means  and  variances 
as  a  means  to  estimating  limits  on  the  s amp le-to-s ample  variations.  It  must  always 
be  kept  in  mind  that  it  is  simply  impossible  to  predict  the  energy  density  spectrum 
of  a  single  transient  with  simultaneous  arbitrary  precision  in  spectral  resolution 
and  amplitude.  This  is  particularly  true  in  cases  where  the  effective  width  of  the 
envelope  spectrum  is  significantly  larger  than  that  of  the  noise  process,  since 
increasing  the  analysis  bandwidth  does  not  produce  a  decrease  in  the  standard 
deviation  of  the  results  until  this  bandwidth  becomes  commensurate  with  the  effect¬ 
ive  spectral  width  of  the  transient  itself. 

In  monopulse  signal  processing  situations  where  either  the  signal  of  interest, 
or  the  interfering  noise,  can  be  modeled  as  a  random  transient  of  the  type  assumed 
in  this  work,  the  results  found  here  predict  the  extent  to  which  filter  output 
energies  vary  from  pulse  to  pulse.  Admittedly,  narrow  band  filtering  represents 
but  one  of  many  processing  techniques  that  might  be  employed.  Still,  it  is  pro¬ 
bably  the  most  basic  and  indicates  the  kind  of  effect  to  be  expected  in  more 
general  situations.  In  predicting  the  signal-to-noise  ratio  at  the  output  of  a 
narrow  band  filter,  or  in  predicting  the  results  of  thresholding  operations 
following  such  filtering,  the  statistical  deviation  avm  sample  to  sample  must  be 
considered.  If  the  system  is  tailored  too  closely  to  the  expected  value  of  the 
output  energies,  it  may  well  malfunction  when  those  samples  arrive,  which,  because 
of  the  inherent  statistical  nature  of  the  ensemble,  yield  output  energies  signi¬ 
ficantly  different  from  the  mean.  Obviously,  this  difficulty  can  be  alleviated  by 
increasing  the  filter  bandwidth  and  thus  decreasing  the  statistical  spread,  but 
often  this  deteriorates  the  signal-to-noise  ratio  out  of  the  device,  and  another 
difficult  nroblem  of  trade-offs  appears.  For  reasonable  statistical  stability, 
the  product  of  the  effective  pulse  duration  and  analyzing  bandwidth  should  be 
appreciably  greater  than  unity,  but  admittedly  this  may  not  always  be  possible. 

It  should  be  noted  that  the  absolute  lower  limit  for  the  output  energy  standard 
deviation  as  a  percentage  of  the  mean  is  on  the  order  of  lA'TW^"  no  matter  what 
filter  bandwidth  is  used.  For  certain  random  transient  ensembles,  therefore, 
particularly  those  whose  effective  spectral  width  is  controlled  by  the  envelope 
spectrum  (and  thus  is  on  the  order  of  1/T) ,  the  output  standard  deviation  will 
always  be  large,  no  matter  how  it  is  analyzed.  Ultimately,  the  use  of  narrow 
band  filtering  must  be  approached  with  care  in  moncpulse  processing  systems  with 
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particular  attention  to  the  opposing  demands  of  finer  resolution  (and  enhanced 
signal-to-noise  ratio)  and  statistical  stability.  This  study  provides  some  in¬ 
sight  and  theoretical  guidance  in  seeking  an  intelligent  compromise. 

Suggestions  for  Future  Research 

There  are  several  avenues  of  future  research  in  both  theoretical  and  experi¬ 
mental  directions  that  can  be  proposed  at  the  present  time.  In  the  latter  area, 
it  would  be  of  great  interest  to  generate  an  experimental  ensemble  of  random 
transients  and  to  subject  representative  sample  functions  to  narrow  band  filtering 
in  an  attempt  to  test  the  predictions  of  the  present  theory.  These  experiments 
could  be  very  efficiently  performed  using  digital  computer  signal  processing 
techniques.  In  particular,  after  producing  a  suitable  random  transient  ensemble 
and  recording  its  members  on  magnetic  tape,  it  would  be  straightforward  to  sample 
and  digitize  the  sample  functions  for  computer  analysis  using  digital  filtering 
and  Fast  Fourier  transform  techniques.  This  would  provide  an  immediate  check  on 
the  appropriateness  of  the  assumptions  and  simplifications  used  throughout  and 
would  lend  a  good  deal  of  insight  to  the  practical  implications  of  sample-to-sample 
variations.  Ultimately,  these  methods  could  be  used  to  study  the  statistical/ 
spectral  properties  of  transient  noise  phenomena  from  the  real  world,  such  as  EER 
returns  and  sonar  reverberation.  This  would  help  in  evaluating  the  model  chosen 
and  in  selecting  appropriate  parameters  for  characterizing  the  signals  of  interest. 

It  should  be  added  that  it  may  well  be  feasible  to  generate  random  transient  sample 
functions  entirely  within  the  computer,  using  digital  filtering  techniques  on  a 
table  of  random  numbers,  and  thus  doing  away  with  any  "analog'1  portions  of  the 
experiment  altogether.  It  remains  to  be  seen  which  approach  will  be  the  most 
effective. 

Many  interesting  questions  have  been  raised  on  the  theoretical  front.  As  the 
reader  has  certainly  noticed,  the  final  results  of  the  study  have  come  out  of  a 
gross  simplification  of  the  spectral  covariance  surface.  The  proposed  experimental 
study  may  well  show  that  this  approximation  is  inadequate,  and  in  that  event,  a 
better  characterization  will  be  necessary.  One  is  faced  with  the  problem  of  devising 
expressions  for  the  form  of  the  surface  which  retain  satisfactory  accuracy  while 
remaining  in  the  realm  of  mathematical  tractability .  Perhaps  an  approach  along  the 
lines  of  an  expansion  in  orthogonal  functions  would  be  appropriate  here. 

The  narrow  band  filter  is  but  one  example  of  a  linear  system,  albeit  one  of 
great  practical  and  theoretical  importance.  It  would  be  desirable  to  extend  the 
present  theory  to  embrace  more  general  linear  systems  with  transient  random  inputs. 

In  particular,  the  correlation  properties  of  the  input  and  output  of  such  systems 
and  the  notion  of  optimum  linear  filtering  when  signal  or  noise  or  both  are  random 
transients  are  of  both  practical  and  academic  interest.  The  application  of  double¬ 
frequency  Fourier  transforms  (see  reference  (h) ,  Chapter  9  and  reference  (m) , 

Chapter  12)  is  another  area  that  bears  study  since  there  is  a  strong  resemblance 
between  that  body  of  techniques  and  those  developed  here. 

The  study  reported  above  has  attempted  to  illuminate  certain  aspects  of  the 
behavior  of  transient  random  signals.  In  particular,  it  has  concentrated  on  the 
problems  that  arise,  due  to  the  pulse  nature  of  the  waveforms,  in  seeking  to 
formulate  a  meaningful  description  of  the  energy  density  spectrum  of  such  signals. 
A^ong  the  most  interestin'?  findings  hnve  v>een  thp  uncertainty  relations  that  exist 
between  spectral  amplitude  and  frequency  resolution  and  th^ir  implications^  in 
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practical  measurement  and  processing  situations.  This  is  but  one  aspect  of  a 
very  general  problem  in  scientific  research;  that  of  defining  and  coping  with  the 
inherent  uncertainties  and  imprecision  of  scientific  measurements.  Certainly 
the  realization  that  there  are  definite  limits  on  our  ability  to  resolve  the 
detailed  properties  of  natural  phenomena  has  been  one  of  the  touch-stones  of 
modern  science,  leading  to  overall  average  statistical  descriptions.  For  this 
reason,  it  is  not  inappropriate  to  close  with  these  perceptive  words  of 
Max  Born  (reference  (j)): 

"The  concept  of  chance  enters  into  the  very  first  steps 
of  scientific  activity  in  virtue  of  the  fact  that  no  observa¬ 
tion  is  absolutely  correct.  I  think  chance  is  a  more 
fundamental  conception  than  causality;  for  whether  in  a 
concrete  case  a  cause-effect  relation  holds  or  not  can  only 
be  judged  by  applying  the  laws  of  chance  to  observations." 


EDWARD  C.  WHITMAN 

Magnetics  and  Electrical  Division 
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APPENDIX  A 

The  Spectral  Normalization  Problem 


As  noted  in  the  main  text,  the  purpose  of  the  spectral  normalization  is  to 
divorce  the  results  of  the  analysis  from  specific  parameter  values,  insofar  as 
this  is  possible.  For  this  reason,  frequency  variables  have  been  expressed  as 
ratios  of  actual  frequency  to  the  "typical"  frequency  of  the  underlying  noise 
process  (u>2  in  the  broad  band  case;  oj^  in  the  narrow),  and  the  length  of  the 

envelope  has  been  expressed  in  terms  of  the  number  of  "typical"  noise  periods 
contained  in  it.  Where  needed,  time  variables  are  normalized  by  division  by  the 
envelope  length  T.  As  an  example,  consider  the  rectangular  pulse  of  broad  band 
noise  for  which 
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From  Parseval’s  Theorem  it  is  known  that 
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Now  introducing  the  parameters 
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and  substituting  these  into  equation  (A-2)  yields  a  function  of  r: 
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Making  a  change  of  variables  in  the  second  integral  yields 
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For  convenience,  we  ask  that  the  normalized  average  total  energy  be  unity. 
Evidently  this  is  attained  for  the  function 


2N  T  E[$'ss(r^ 
o 


(A-S 


which  gives 


¥*3s<r>:>  - 


(0  T 
2 


(1-x) 


sin2IIqx 

2Jlqx 


cos  2IIqxr  dx 


(A- 1C 


or  since  q  = 


v 

2n 


E»t*  (r)] 


2q  1  ^-X>  C0S  2IIqXr  ^ 


(A-l] 
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It  is  easily  seen  from  above  that 


00 

f  EN[*ss(r)]  dr  =  1 

—  00 


and  that 


E|>gg(w)] 


(A-12) 


(A-13) 


The  extension  of  this  normalization  to  the  covariance  expressions  is  straight¬ 
forward  and  results  in  a  definition  of  CovN[$  (r),  <J>  (s)]  such  that 

Cov[$ss(u,)?  *ss((o‘)]  =  4N02T2  CovN[^s  (  =  $ss(s  -S7)] 

One  goes  through  a  similar  procedure  to  obtain  the  analogous  expressions  for  the 
narrow  band  cases. 

There  is  one  disadvantage  to  the  present  normalization  which  may  mislead  the 
unwary  unless  it  is  pointed  out.  We  know  that  E[$gs(w)]  is  the  Fourier  transform 

°f  E[4*ss(t)]  and  vice-versa.  What  can  be  said  of  the  Fourier  transform  of 

Ew[$  (r)]?  From  equation  (A-ll)  one  can  see  that 
in  s  s 


=  q  f  u-M) 


-1 


2nq 

-  Iff  ;  (1 


sin2llqx  -j2;[qrx 
2Hqx  6 


dx 


(A-15) 


x  v  sin  x  -jxr 


2IIqV  x 


dx 


-211  q 


Thus  EM[$  _(r)]  and  the  function 

W  SS 


«  -  C1  -  M)  ^ ,  -2n<i  - x  - 2,ui 


(A-16) 


0  ,  elsewhere 


A-3 
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are  a  Fourier  transform  pair.  This  is  not  the  expression  one  gets  from  equation 
(A-l)  by  setting  2N0u2T  =  1  and  substituting  the  parameters  of  equations  (A-4), 
As  a  matter  of  fact,  under  these  conditions,  equation  (A-l)  becomes 


E[*ss(x)J 


(1  -  1*1) 


sin2IIqx 

2IIqx 


1  <  x  <  1 


(A-17) 


0  ,  elsewhere 


and  thus  is  not  the  Fourier  transform  of  (r)  ] . 
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APPENDIX  B 

THE  "WIDTH"  OF  THE  CONVOLUTION  OF  TWO  FUNCTIONS 


In  the  main  body  is  made  the  statement  that  in  an  approximate  sense,  "the 
width  of  the  convolution  of  two  even  functions  is  equal  to  the  sum  of  the  widths 
of  the  functions  themselves."  This  approximation  sees  constant  use  in  practical 
treatments  of  modulation,  transient  analysis,  and  sampling  theory.  This  appendix 
gives  explicit  justification  for  the  cases  of  major  interest  in  this  report. 

If  f(x)  and  g(x)  are  two  real  functions,  their  convolution,  subject  to  certain 
existence  conditions,  is  given  by 


c(x)  f(x)@g(x)  =  f  f(u)  g(x-u)du  =  /  g (u)  f(x-u)du 


(B-l) 


If  f(x)  and  g (x)  are  both  even  functions,  then 


g(x-u)  =  g(u-x) 


(B-2a) 


and 


f(x-u)  =  f(u-x) 


and  thus 


(B-2b) 


00  00 

c  (x)  “  /  f  (u)  g(u-x)du  =  /  g(u)  f  (u-x)du  =  <J>,  (-x) 

-00  _oo 


(B-3) 


where  <J>fg(x)  is  the  cross-correlation  function  of  f(x)  and  g(x)  defined  as 

CO 

$£g(x)  =  f  f(u)  g(u+x)du  (b-4) 

—00 


Furthermore,  when  both  functions  are  even, 


^fg(x)  =  *fg(~x) 


and  thus 


c(x)  -  4>fg(x) 


(B— 5  ) 


B-l* 
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which  is  to  say  that  the  convolution  of  f(x)  and  g(x)  is  the  same  as  their  cross¬ 
correlation  function.  In  the  present  application,  we  deal  with  power  and  energy 
density  spectra,  which  are  indeed  real,  even  functions,  integrable  over  the  whole 
real  line.  In  the  low  pass  cases,  these  functions  are  concentrated  at  the  origin 
and  tend  to  zero  as  their  arguments  approach  +«>  or  -°°  .  One  can  thus  define  a 
spectral  half  width  W  which  represents  the  positive  frequency  range  for  which  the 
spectral  amplitude  differs  significantly  from  zero.  In  band  pass  situations,  a 
spectral  width  can  similarly  be  defined  for  the  spectral  concentration  around  the 
center  frequency.  Admittedly,  the  criteria  one  uses  to  define  the  spectral  width 
are  arbitrary  and  subject  to  question,  but  if  it  is  consistently  defined  in  the 
problem  of  interest,  the  interpretation  should  be  unambiguous. 

As  an  example,  consider  the  case  where  both  convolved  functions  are  rectangular 
pulses  of  unit  amplitude,  entered  at  the  origin: 


f(x) 

=  1, 

-A  £  x  _<  A 

(B— 6 ) 

=  0, 

elsewhere 

g(x) 

-  1, 

-B  _<  x  £  3 

(B— 7) 

=  o, 

elsewhere 

Assuming  that  A  _>  B  and  using  equation  (B-3) ,  we  find  that 

c(x)  ■  f  (x)  @  g  (x)  ■  x  +  A  +  B ,  -  A-  B_<x<-A  +  B 


*  2B  ,  -  A+  B£x£A-B  (B-8) 


=  A  +  B-x,  A-B£x£A+B 

=  0  ,  |  x  |  >  A  +  B 

This  is  shown  graphically  in  Figure  B-l.  If  the  function  half -widths  are  defined 
as  the  positive  frequency  spread  for  which  they  are  non-zero,  it  is  evident  that 


Mj 

It 

> 

(B-9a) 

W  =  B 
g 

(B-9b) 

and  that 


W„  -  A  +  B  =  V*  +  V 
c  f  g 


(B-10) 


B- 


O 
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which  confirms  the  result.  This  is  a  very  important  special  case,  since  in 
practice  one  often  approximates  spectra  of  arbitrary  shape  as  rectangular  idealiza¬ 
tions  whose  amplitude  is  set  equal  to  the  maximum  value  of  the  spectrum  of  interest 
and  whose  "equivalent  width"  is  chosen  to  set  the  area  under  the  approximation  equal 
to  that  under  the  actual  curve. 


Another  example  of  interest  is 

that  where  f(x)  and  g(x) 

are  both  double 

sided  decaying  exponentials: 

f(x>  =  e-H/A 

8 

V 

X 

V 

8 

1 

(B-ll) 

/  \  -jx|/B 

g(x)  =  e 

8 

V 

X 

V 

8 

1 

(B-12) 

with  A  >  B.  A  reasonable  measure  of  the  function  half -width  in  this  case  is  that 
point  where  the  function  amplitude  falls  to  1/e  times  its  value  at  the  origin. 
This  gives 

Wf  =  A  (B-13a) 

Wf>  =  B  (B-13b) 

From  equation  (B— 3) ,  we  find  that 


c (x)  =  f (x)  (^g(x) 


AB  (  -M/a 

A  +  B  \e 


+ 


AB 

A  -  B 


(y|x|/A 


-Ul/B) 

e  '  > (B-13) 


—oo  <  X  <  00 


The  constituent  functions  and  the  resulting  convolution  are  shown  for  two  cases  in 
Figure  B-2.  If  the  half-width  of  the  convolution  is  taken  as  the  argument  for 
which  the  amplitude  falls  to  1/e  times  its  value  at  the  origin,  the  following  table 
of  computed  results  can  be  derived  using  equation  (B-13) : 


t-h 

II 

> 

W  =  B 

R 

Wr  +  W 
f  R 

W  =  c  1  (1/e) 
c 

1.0 

0.9 

1.9 

2.06 

1.0 

0.5 

1.5 

1.59 

1.0 

0.25 

1.25 

1.29 

1.0 

0.1 

1.1 

1.11 

B-3 
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One  again  concludes  that  Wc  & Wj  +  Wg  with  the  approximation  becoming  better  and 
better  as  g(x)  becomes  increasingly  narrow  with  respect  to  f(x). 

The  generality  of  this  conclusion  is  made  even  more  plausible  by  those  inter¬ 
pretations  of  the  convolution  integral  which  treat  g(-s)  as  a  "scanning  function" 
which  slides  along  f(x)  to  generate  a  smoothed  version  of  the  latter  (see  refer¬ 
ence  (k) ,  p  323,  and  reference  (1),  p  71).  Naturally,  the  wider  g(x)  becomes,  the 
more  broadly  will  f(x)  be  smoothed  out,  and  the  width  of  the  smoothed  version 
will  be  proportional  to  the  sum  of  the  widths  of  f(x)  and  g(x). 

It  is  easily  seen  in  this  graphical  interpretation  that  in  cases  where  one  of 
the  two  convolved  functions  is  of  band  pass  nature  and  the  other  of  low  pass 
character,  the  resulting  convolution  will  yield  a  smoothed  (broadened)  version  of 
the  former  centered  at  the  original  center  frequency  with  an  equivalent  width 
approximately  equal  to  the  sum  of  the  double  sided  width  of  the  two  constituents. 
This  is  precisely  what  is  done  in  the  familiar  spectral  interpretation  of  amplitude 
modulation. 


B-4 
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i  i  i  i  i  i 

-A-B  -A+B  1  A-B  A+B 

FIG.B-1  THE  CONVOLUTION  OF  TWO  RECTANGULAR  PULSES. 


B-5 


FIG.B-2  CONVOLUTION  OF  DOUBLE-SIDED  DECAYING  EXPONENTIAL  PULSES  FOR  TWO  SETS 
OF  PARAMETER  VALUES.  ALL  AMPLITUDES  HAVE  BEEN  NORMALIZED  TO  UNITY  AT 
THE  ORIGIN. 
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APPENDIX  C 


An  Alternative  Approach  to  the  Narrow  Band  Filtering  Problem 
with  a  Solution  for  an  Important  Special  Case 


Consider  passing  a  random  transient  ensemble  through  the  narrow  band  filter 
given  by  equation  (109)  and  Figure  31.  By  equation  (113), 


w^B/2 

E[F(H)  ]  =  2  /  E[  *__(<■>)]  du  (C-l) 

o)c-B/2 


Now  by  equation  (8)  this  is 


E[F(H)] 


n 


u)c+B/2 


»c-B/2 


/  ljl  (x)  (j> 
rnn  Y 


ee 


-jwx 

(x)e  dx  dm 


( C— 2) 


Interchanging  the  order  of  integration  yields 


E[F(H) ] 


oo  toc+B/2 

_  _JL  /  ij;  (x)  <j>  (x)  f  e  ^u3£  dw 

n  -oo  nn  ee  u  -b/2 

C 


dx 


(C— 3) 


_B_ 

n 


.  xB 

00  sm  -r~ 

f  *nn(x)  *ee(x) - ^ -  e"jlocX  dx 


xB 

2 


since  ^nn(x)  and  4>ee(x)  are  both  even  functions  of  x. 


From  equation  (115) ,  we  know  that 


o>c+B/2  b)c+B/2 

Var[F(H)]  =  4  /  /  Cov[4i  s(to),  <*>  (to’)]  do)  dw'  (C-4) 

a)c-B/2  u>c-B/2 
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where  Cov[$gS(w) ,  4>ss (<*>’)]  is  given  by  equations  (41)  and  (44)  as 


Cov[$ss(io),  $S8(a)’)] 


(C-5) 


00  00  00  00 

-ftp  f  f  f  f  e(t)e(x)e(u)e(v)^nn(t-u)i|<nn(x-v)e-Jtl)(t^x)e->,^u-v)dt  dx  du  dv 


0  0  0  0 


+  ~n?z  f  f  f  f  e(t)e(x)e(u)e(v) ik._ ( t-v)  tJl.  (x-u) e-^ u ^ t-x^ e~^u '  dt  dx  du  dv 


4nT 


*'nn"-  '  ynn 


oooo 


When  equation  (C-5)  is  used  in  equation  (C-4) ,  some  judicious  manipulation  makes 
it  possible  to  perform  the  integrations  on  w  and  oj*.  The  following  expression 
emerges : 


00  00  00  00 

VarlF(H)]  =  :p-  f  f  f  f  e(t)e(t-a)e(u)e(u-b) 

r-oo  ^00  0  0 


(C-6) 


^nn^t-u^nn^t_a_u+b^  +  'J'nn^-^Hnn^-11-3^ 


sin  aB/2  sin  bB/2 
aB/2  bB/2 


e  ^wo^a  b^dt  du  da  db 


Now  let  us  consider  the  special  case  where  the  effective  noise  bandwidth  is  much 
larger  than  the  analyzing  filter  bandwidth,  which  in  turn  is  much  larger  than  the 
effective  bandwidth  of  the  envelope  function.  This  implies  that  tk,_(x)  is  very 

much  nattower  than  - B^/2 —  *  wb*cb  *n  turn  is  much  narrower  than  the  envelope 

autocorrelation  functions.  An  approximate  solution  for  this  case  is  found  by 
assuming  that  the  noise  correlation  function  can  be  well  represented  by  the  Dirac 
delta  function: 


(C— 7) 


Then,  from  equation  (B-3) , 


C-2 


NOLTP.  68-124 


E[F(H)]  a? 


~  J  *WX>  *ee<x>  -""'btI'"-  e"jt°cx 


(C-8) 


E[F(H)]  «  iee(0)  =  /  e2(t)  dt 

J1  o 

Substituting  equation  (C— 7)  into  equation  (C-6)  yields: 


N02B2 


CO  CO  CO  CO 


Var [F(H) \tz.  — jp-  f  /  /  /  e(t)e(t-a)e(u)e(u-b)yQ(t-u)y0(t-a-u+b) 


—00  —00  0  0 


sin  aB/2  sin  bB/2  -jo>0(a-b) 


dt  du  da  db 


B2N  2  00  00  “  ” 

+  - 2—  f  f  f  f  e(t)e(t-a)e(u)e(u-b)y  (t-u+b)u  (t-u-a) 

TT  £  u  Q 


—co  —co  o  0 


-~gy~  sl"B)°/2  e-Ma-V  dt  du  da  db 


The  integrand  of  the  first  integral  is  zero  except  where 


(C-9) 


t-u  =  t-a-u  +  b  =  0 


(ClOa) 


and  that  of  the  second  is  zero  except  where 


t-u  +  b  =  t-  u-  a  =  0 


(ClOb) 


This  implies  in  both  cases  that  a  =  b.  Thus,  integrating  on  u  and  a  yields 


Var[F(H)]  S£ 


2M  2B2  0°  oo 


/  /  e2(t)e2(t-b)  -^sLMZ|  dt  db 

-oo  o  (bB/2)2 


(C-ll) 
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where  we  have  used  the  integral  properties  of  the  Dirac  delta  function.  Now  define 
the  autocorrelation  function  of  the  square  of  the  envelope  function: 


00 

<j>  2  2  (T)  =  f  e2(t)e2(t— r)  dt 

e  e 

o 


(C-12) 


and  with  this  at  hand, 


Var[F(H)] 


4N02B2 


/  ^e2e2(b) 

o 


sin2  bB/2 
(bB/2)z 


db 


(C-13) 


Now  recalling  that  _sin  bB/2 _  has  been  assumed  to  be  much  narrower  than  the 

bB/2 

envelope  autocorrelation  function. 


Var[F(H)]  « 


4N02B2 


*e2e2(°) 


sin2  bB/2 
(bB/2)2  ' 


db 


(C-14) 


Consulting  a  table  of  integrals  then  yields 


Var[F(H)] 


4N02B 

n 


<i>e2e2(0) 


4Nd2B 

n 


00 

/  e4(t)  dt 
o 


(C-15) 


Finally,  forming  the  ratio 


R(H) 


/Var  F(H) 
E[F(H)] 


(C— 16) 


which  is  the  normalized  standard  derivation  of  the  output  energy  for  this  case, 
we  find  that 


R(H)  = 


2/ar 

sr 


/ar 


/f  e4(t)dt 


f  e2(t)dt 


( C— 17) 
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When  e(t)  is  a  rectangular  pulse  of  length  T,  the  ratio  of  integrals  becomes 
merely  l/^T~ .  Thus  for  this  case, 


R(H) 


2/T 


(C-18) 


which  is  exactly  the  expression  found  using  the  method  of  spectral  covariance 
surfaces  (see  equation  (143))  and  which  is  quoted  in  reference  (h)  as  arising 
in  a  "well-resolved"  spectrum  analysis.  Of  course,  all  of  the  above  derivation 
after  equation  (C-6)  is  heuristic  rather  than  rigorous,  but  it  provides  an 
interesting  confirmation  of  a  previous  result  and  demonstrates  the  importance 
of  considering  the  relative  size  of  the  several  bandwidths  involved. 
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A  class  of  transient  random  signals  is  modeled  as  the  product  of  a 
deterministic,  square  integrable  envelope  function  and  a  Gaussian  random  process 
having  a  well-defined  power  spectrum.  The  passage  of  an  ensemble  of  such  random 
signals  through  a  linear  filter  is  studied  with  particular  emphasis  on  the  mean 
and  variance  of  the  total  output  energy.  It  is  found  that  an  important  role  is 
played  in  these  considerations  by  the  covariance  function  between  values  of  the 
energy  density  spectrum  of  a  sample  function  evaluated  at  different  frequency 
arguments.  Accordingly,  the  form  of  this  function  is  derived  and  portrayed  as  a 
surface  lying  above  a  two-dimensional  frequency  plane.  Examples  of  these  spectral 
covariance  surfaces  are  presented  and  discussed  for  both  rectangular  and  decaying 
exponential  pulses  of  both  broad  and  narrow  band  Gaussian  noise,  and  their  general 
characteristics  are  identified.  Finally,  the  problem  of  idealized  narrow  band 
filtering  is  specifically  approached  and  approximate  expressions  derived  for  the 
mean,  variance,  and  normalized  standard  deviation  of  the  output  energy  of  a  narrow 
band  filter  excited  by  rectangular  pulses  of  narrow  band  Gaussian  noise.  The 
approximate  relationships  of  the  filter  bandwidth,  pulse  duration,  and  underlying 
noise  spectrum  are  explored  for  their  effects  on  spectral  resolution  and  statistical 
stability,  leading  to  an  uncertainty  principle.  The  implications  of  these  findings 
for  spectral  analysis  and  monopulse  signal  processing  are  discussed  in  the  light  of 
this  uncertainty  principle  and  the  limitations  it  imposes  on  the  simultaneous 
precision  of  frequency  resolution  and  spectral  amplitude. 
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